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ABSTRACT 

The subject of relativistic hydrodynamics is explored using the tools of gauge/gravity 
duality. A brief literature review of AdS/CFT and gauge/gravity duality is presented 
first. This is followed by a pedagogical introduction to the use of these methods in 
determining hydrodynamic dispersion relations, w(q), of perturbations in a strongly 
coupled fluid. 

Shear and sound mode perturbations are examined in a special class of gravity duals: 
those where the matter supporting the metric is scalar in nature. Analytical solutions 
(to order and respectively) for the shear and sound mode dispersion relations are 
presented for a subset of these backgrounds. 

The work presented here is based on previous publications by the same author |I] , [2] , 
and though some previously unpublished results are also included. In particular, the 
subleading term in the shear mode dispersion relation is analyzed using the AdS/CFT 
correspondence without any reference to the black hole membrane paradigm. 
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of the free parameter C2 appearing in the metric (4.76). Each plot corre- 
sponds to a different number of spatial dimensions p. This range of values 
of C2 was chosen because outside this particular range, unphysical behav- 
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Chapter 1 



Introduction 



Each of the fundamental forces in nature is currently understood in the context of 
a particular theory. Gravity is understood in terms of Einstein's theory of general 
relativity (the quantum theory of gravity is currently unknown, and highly sought after); 
while the electromagnetic and weak forces are both understood in the context of the 
standard model of electroweak interactions. The subject of this research is the theory 
of the strong nuclear force, Quantum Chromodynamics (QCD). 

QCD is the theory of the interaction that binds protons and neutrons together to 
form atomic nuclei. Quarks (fermions), and gluons (gauge bosons of the underlying 
SU(3) gauge group) are the fundamental constituents of the theory. It is believed that 
if one heats normal matter to high enough temperatures, or compresses it to high enough 
densities, a new state of matter is formed wherein the boundaries between particular 
nucleons are no longer well defined. Instead of a sea of nucleons, one finds a "soup" of 
quarks and gluons. We will refer to this state of matter as quark-gluon plasma (QGP) 
It is believed to have dominated the universe in the first few microseconds after the Big 
Bang. 

Experimentalists are attempting to re-create this environment using particle accel- 
erators which smash nuclei into each other at very high energies. In the last few years, 
there is evidence that the Relativistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Laboratory has created such a state of matter IH [5l |6l [7] . 

The plasma created at RHIC exhibits some surprising properties. In particular, 
observation of strong collective behavior (elliptic flow), and the large energy loss of 
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high energy particles traversing the medium (jet quenching) indicate that the plasma 
interacts very strongly with itself, and is thus referred to as strongly coupled. This 
presents a problem for theoretical physicists, because the equations of QCD cannot 
be solved analytically in this regime; traditional approximation techniques involve a 
perturbative expansion in the coupling constant which measures the strength of the 
interaction. If the coupling constant is large (i.e. 0(1)), such a perturbative expansion 
is unreliable, as each successive term in the series becomes larger than the previous 
one, and the approximation scheme breaks down. This difficulty with QCD has plagued 
theorists for over 30 years since the inception of the theory. 

However, in the late 1990's, new techniques were developed which attempt to address 
these non-perturbative problems [H |9l [10] . In particular, it has been conjectured that 
there is a duality between certain strongly coupled gauge theories and weakly coupled 
string theories. By this, it is meant that both theories may describe the same physics, 
but calculations may be easier in one theory than the other. As mentioned above, it 
has been argued that the duality is a weak/strong type, such that when one theory 
is strongly coupled, the other is weakly coupled and vice-versa. One can immediately 
see the possible application to QCD; if a dual theory to QCD is found, one can do 
computations in the dual theory (at weak coupling, where calculations are possible), 
instead of in QCD itself (at strong coupling where calculations are difficult). 

Much work needs to be done before we can apply these ideas to real- world QCD. 
At the present time, the duality (referred to in the title as the gauge/gravity correspon- 
dence) is well developed for a theory which shares some properties with QCD, but lacks 
some of the essential features that we observe in the real world. Attempts are currently 
being made to both search for a dual theory of QCD, and to modify the idealized theory 
to make it more physically relevant. 

The ultimate goal of this line of research is the calculation of observable properties 
of the quark-gluon plasma. This is still a difficult task, because the string dual to 
QCD is not currently known. In light of this fact, there are two avenues which one may 
pursue. First, one can attempt to compute observables in phenomenologically motivated 
extra-dimensional models of QCD. The drawback here is that at the present time, such 
models are lacking from a theoretical point of view. Secondly, one can attempt to look 
for universal features of strongly coupled theories, which might then also be applicable 
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to QCD. 

Here we employ the latter approach; we search for universal behavior in hydrody- 
namic transport coefficients which describe a fluid's response to small perturbations. 
Hydrodynamics is an effective theory of fluids in thermal equilibrium. It is applicable 
when macroscopic length and time scales of interest are much longer than any micro- 
scopic scale. In this regime, the effective stress-energy tensor can be constructed as 
a derivative expansion in the fluid velocity. Naturally, such an effective construction 
introduces unknown coefficients. These are the transport coefficients mentioned above, 
and which will be discussed in much greater detail in the coming chapters. 

These coefficients are necessary input for hydrodynamical simulations of the quark- 
gluon plasma. Currently, we understand that a heavy ion collision occurs in roughly 
three distinct phases. Before the collision, the two colliding nuclei are highly Lorentz 
contracted, but have not yet entered into thermal equilibrium; in this phase, the nuclei 
can be described within the framework of, for example, the Color Glass Condensate [11] . 
After the collision, some of the energy of the colliding nuclei is transformed into a bath 
of new particles, which rapidly thermalizes. Once the new system has reached thermal 
equilibrium, a hydrodynamic description is appropriate. The newly created region of 
quark-gluon plasma then expands and cools; once it has cooled enough, it drops out of 
thermal equilibrium and the matter hadronizes into particles which are then detected. 
Describing the initial (pre-equilibrium) and final (freeze out) stages of the heavy ion 
collision is a difficult and very active area of research, it is also beyond the scope of 
this thesis - for reviews, see [121 [131 E]- Luckily, it appears that there is a stage of 
the heavy ion collision when the matter is thermally equilibrated, and thus admits a 
hydrodynamic description. It is in this phase of the collision that one needs to specify 
the relevant transport coefficients. It is thus desirable to calculate these coefficients, but 
as mentioned above, the strong coupling of the plasma renders traditional perturbative 
calculations unreliable. Thus, one is led to try using gauge/gravity duality methods to 
compute such quantities. 

The thesis is organized as follows 

• In Chapter [2j we give a brief review of the motivation, and some of the basic tools 
of gauge/gravity duality. While hydrodynamics and transport coefficients are the 
main focus of this work, in this section we also discuss some other applications of 
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gauge/gravity duality to strongly coupled plasma. Particular attention is paid to 
thermodynamic applications, some of which are useful in later chapters. 

• In Chapter [3j we explain the relevant background material regarding hydrody- 
namics. We also discuss some of the standard methods of computing transport 
coefficients from gauge/gravity duality, and explain in detail the method which 
will be used in deriving the central results of the thesis. 

• In Chapters |4] and [5} the techniques described above are applied to a particular 
class of dual theories. These sections contain the central results of this thesis. We 
derive equations which are applicable to a wide variety of gravitational dual theo- 
ries, and could in principle be used within the context of phenomenological models 
of QCD. We also apply these equations to a few analytically solvable special cases, 
and examine the results to determine if any universal features of strongly coupled 
plasmas are manifest. The methodology employed in each of these chapters is the 
same, but each one concerns a different type of hydrodynamic perturbation. In 
Chapter |4] we study shear perturbations while sound (or compressional) pertur- 
bations are studied in Chapter [5j 

• Finally, in Chapter |6] we discuss the main conclusions of the analysis and also 
mention open questions and prospects for future investigation. 



Chapter 2 

Gauge/Gravity Duality: An 
Overview 

In this chapter, we first explain the basics of gauge/gravity duality, and the AdS/CFT 
(Anti-de Sitter /Conformal Field Theory) correspondence. The focus of this thesis is 
phenomenology, and particular applications of the correspondence, so many of the de- 
tails regarding the foundations of the duality will not be presented here. Relevant 
references are provided for the interested reader. Some excellent reviews of the subject 
which parallel the present discussion can be found in \15\ [TBI [T71 \T8\ [T9] . 

Once the basics of the duality have been established, we present a review of the 
relevant literature regarding applications of the duality (mass spectra, form factors, jet 
quenching, etc.). The main application discussed in this thesis (hydrodynamic disper- 
sion relations and transport coefficients) will be explored in more detail in the next 
chapter. 

Finally, we discuss in some detail applications of the duality to thermodynamic 
aspects of the relevant gauge theory. Many of these results will be useful in subsequent 
chapters. 

2.1 What is Gauge/Gravity Duality? 

The original motivation for gauge/gravity duality came from considerations of type IIB 
string theory, which exists in ten dimensions. It is not necessary to understand the 
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details of string theory for this derivation, but one should appreciate that it is a theory 
which contains strings, but also extended membrane like objects, branes. D-branes are 
such objects which can serve as endpoints for the strings ('D' stands for Dirichlet). The 
theory has two types of excitations, closed strings which have no endpoints, and open 
strings which are allowed to have endpoints on a D-brane. Strings are characterized by 
their length Ig, and the coupling constant which controls the strength of their interaction 
is denoted gs- The string coupling can be related to the D dimensional Newton's 
gravitational constant Gd by 

Go^g'sl?-'. (2.1) 
To proceed, let us consider type IIB string theory in the background of a system of 
N D3-branes stacked on top of one another. (A D3-brane is a D-brane which extends in 
three spatial dimensions). We will now consider the low energy dynamics of this system 
from two different viewpoints. 

First, the low energy excitations of open strings on the system of N D3-branes can 
be described by a SU{N) gauge theory. Specifically, the relevant low energy effective 
theory is AA = 4 supersymmetric Yang-Mills (hereafter SYM) theoryj^ The gauge 
coupling qy m is related to the string coupling gg as [20] , 

g^M = ^T^9s- (2.2) 

If we consider the string length to be small, then we can neglect any interactions between 
the branes and the closed strings which exist outside of the branes (in the region called 
the hulk). We can also neglect any interactions among the closed strings themselves. 
This is because all interaction terms will contain positive powers of the gravitational 
coupling Gd. If we are interested in processes which have typical energy E much smaller 
than the Planck mass Mp, these interaction terms will be suppressed by powers of E /Mp. 
Thus, from this point of view, the low energy limit of this entire system consists of two 
decoupled pieces: free low energy closed string theory (type IIB supergravity) in the 
bulk, and = 4 SYM theory on the branes. This point of view is presented in graphical 



form in Fig. 2.1 



Now, consider this system from an alternative point of view. If N is very large, the 
stack of branes will have a considerable amount of energy, and thus will curve space-time 



^ One should take care to distinguish J\f (the number of supercharges in the supersymmetric gauge 
theory) from A'^ (the number of D-branes and the number of colors in the gauge theory). 



Viewpoint 1 
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Figure 2.1: Type IIB string theory in the background of a set of N coincident D3 branes. 
Low energy excitations of this system decouple into the two pieces shown above; the 
low energy excitations of the brane system are those of a gauge theory. 



in accordance with general relativity. In [21 , it was shown that one can find black hole 
type solutions to these supergravity equations. The classical metric which solves the 
supergravity equations can be written 



1 



L4 



dt"^ + dx^) + V 1 + ^ (dr^ + r'^d^l) , 



(2.3) 



where L is the only scale, and is the so called curvature radius. The vector x runs over 
the three spatial coordinates, and 0,^ is the five dimensional angular element. In addition 
to this metric, there is a scalar field (dilaton), and a five- form, but these additional fields 
are unimportant for the argument here. From this point of view, the system now has 
two types of low energy excitations: closed string excitations which are far away from 
the brane which is located at r = 0, or any sort of excitation near the location r = 0. 
From the point of view of an observer at r = oo, these excitations appear to be low 



energy due to the gravitational redshift from the metric. In the near horizon regime 
r « L, the metric becomes 

ds^ = ^ {-dt^ + dx^) + ^dr^ + L^dnl (2.4) 

which is the product of five dimensional anti-de Sitter space, with a 5-sphere {AdS^ x 55). 
(Anti-de Sitter space is a maximally symmetric metric which is a solution to Einstein's 
equations with a negative cosmological constant). 

Furthermore, both of the types of excitations listed above decouple if the energy is 
low enough, because the wavelength of the closed string excitations will be larger than 
the gravitational size of the branes (~ L). This point of view is presented graphically 



in Fig. 2.2 



Now let us compare the two viewpoints of this system. In each description there are 
two decoupled pieces. In the first system, the decoupled pieces are the = 4 SYM gauge 
theory, and low energy closed string excitations. In the second description, the decoupled 
pieces are the excitations in the near horizon regime of the black hole, and low energy 
closed string excitations far away from the horizon. Maldacena's conjecture is that these 
two systems should describe the same physics |8! . Each of the two descriptions has low 
energy closed string excitations as one element, thus we should equate the remaining 
two elements. The result is the conjecture that AA = 4 SYM theory is dual to type 
IIB string theory on AdSs x S5. 

To see why this is useful, consider describing the type IIB string theory by a classical 



supergravity metric (2.4) and associated supergravity fields. In other words, we neglect 
all stringy effects and treat the theory classically. When is such an approximation valid? 
Clearly, one must have the string length much less than the curvature radius: 

< L. (2.5) 

We can relate this to the parameters of the gauge theory by examining the total en- 
ergy/volume of the stack of branes within these two different viewpoints. 

On one hand, the total energy/volume of the brane system is just N times the 
tension of one brane. The branes extend in three spatial dimensions, thus by dimensional 
analysis, 

M/V ~ A//^. (2.6) 
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Figure 2.2: The same system as in Fig. 2.1 viewed in a different way. If the stack of 
branes is large, they curve space-time and can be described by a classical metric. The 
gravitational redshift causes excitations near the horizon to appear very low in energy. 



Clearly, the tension must scale linearly with N ^ and the string length is the only relevant 
scale available to get the dimensions correct. 

On the supergravity side, one can compute the total mass of the D-brane system by 
using the ADM mass formula 

M~y" d^xT™, (2.7) 

where D is the number of dimensions, and T^^ is a particular component of an energy- 
momentum tensor. (We are not being rigorous here, we are doing dimensional analysis 
only). Now, by Einstein's equations, we have 



(2.8) 
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where G^i^ is the Einstein tensor. The components of the Einstein tensor in question 
win in general be functions of the space-time coordinates, but the ADM mass formula 
integrates over these coordinates. The only scale that can enter is the curvature scale 
L. Recalling that in D dimensions, Newton's constant has mass dimension D — 2, by 
dimensional analysis 

M ~ L^-^Gd, (2.9) 
M/V ~ L^-^/Gd- (2.10) 

Furthermore, Newton's constant can be related to the string length, since the super- 
gravity description is only a low energy limit of string theory. By dimension, 

Gd ~ /f-', (2.11) 

finally, 

M/V ~ L^-^/l^-^. (2.12) 



Specifying to D = 10, and equating this result with (2.6), we have 

Nr^{L/ls)^. (2.13) 

The above analysis was done on dimensional grounds only, which makes it impossible 
to track dimensionless quantities like the string coupling Qs . After a more careful analysis 
|22l [23| one finds that the brane tension goes as gj^, and that the gravitational 
constant Gd goes as as mentioned previously. Including these factors, we have the 
useful relation 

gsN ~ g^j^jN ~ L^/li (2.14) 



Here we have used (2.2) in the first step. Going back now to (2.5), one finds that the 
supergravity description is valid for Ig <^ L which in turn implies 

qImN » 1. (2.15) 

The quantity gyj^^jN is referred to as the t'Hooft coupling. We have just shown that 
our classical approximation is valid when dual gauge theory is strongly coupled! 

We also need to suppress quantum gravitational corrections in order to use the 
classical metric. To do so, we require that 

L > (2.16) 
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where Ip is the Planck length. It is related to Newton's constant as 



Gd ~ Ip 



D-2 



(2.17) 



Thus we have the relation 



(2.18) 



Hence, we can neglect quantum corrections provided that 



7V» 1. 



(2.19) 



This is referred to the t'Hooft limit. It is the limit where the rank of the gauge group is 
large N ^ 1, and the t'Hooft coupling is fixed and large Sym^ This is the limit 

where the supergravity approximation is valid. 

To summarize, Maldacena's conjecture is that a certain kind of string theory on 
a curved AdS background is equivalent to A/" = 4 SYM theory. If we neglect stringy 
effects, we can describe the string theory by a classical metric (and other classical fields) . 
It turns out that the regime when this works is the regime when the gauge theory is 
strongly coupled. Hence, we see the usefulness of this approach; the strong coupling 
regime of the gauge theory is exactly the regime where it is difficult to use traditional 
perturbative techniques. However, this is precisely the regime in the dual string theory 
where calculations are easy, because the theory can be described classically! 

We will discuss how one can use the dual gravity theory to do calculations in future 
sections. In the next section, we will discuss how such a duality may be relevant to 



2.2 Quantum Chromodynamics and AdS/CFT 

A phenomenologist would like to apply these techniques to QCD, a real theory of nature, 
as opposed to the supersymmetric Yang-Mills theory mentioned in the previous chapter. 
Are these theories so different? 

For one thing, A/^ = 4 SYM theory is a conformal field theory. Such a theory has 
no intrinsic energy scale (unlike QCD, where the running of the coupling introduces 
the scale Aqcd)- Thus the SYM theory has a gauge coupling which does not run, and 
the theory is not confining. It is also supersymmetric, and has AT S> 1 whereas QCD 



QCD. 
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has = 3. Furthermore, QCD contains quarks whereas the only fermions in the SYM 
theory are the super-partners of the gauge bosons. These fermionic partners transform 
in the same way as the gauge bosons, but this is not how quarks transform. (In more 
pedestrian language, the gauge bosons have two "color indices", as opposed to quarks 
which only have one). 

At first sight, these two theories are thus very different, and one may be skeptical 
about the application of this duality to QCD. In general, there are three philosophies 
regarding the application of gauge/gravity duality to QCD. 

2.2.1 Top Down Approach 

The first approach, usually referred to as the top down approach, focuses on attempting 
to discover other dualities where the gauge theory more closely resembles QCD. Here, 
one must modify the existing string theory in some way. For example, one can address 
the lack of quarks by introducing 'flavor branes' in the bulk |25| [26] . One must break 
conformal symmetry in order to introduce confinement; some examples of string duals 
which do so are those of Witten j27|, Polchinski and Strassler |28j, Maldacena and 
Nunez [29] , and Klebanov and Strassler [30] . The addition of flavor into Witten's model 
by Kruczenski et al. [31^ and by Sakai and Sugimoto |32j allows for chiral symmetry 
breaking. These latter constructions share many features with QCD and are sometimes 
called Holographic QCD. 

2.2.2 Bottom Up Approach 

The second approach (sometimes called the bottom up approach, or AdS/QCD) exam- 
ines the duality from a phenomenological perspective. The idea is to modify by hand 
the existing classical AdS^ background in order to introduce some essential features of 
QCD. The resulting modified backgrounds may or may not be embeddable into string 
theory. The benefits of such an approach are twofold. First, one may view the resulting 
background as a model of QCD, which may be systematically improved and can be 
used to make predictions, even though this model was not derived from a fundamental 
theory. Second, by examining the properties of the backgrounds which lead to QCD 
like behavior, one might gain information that may be useful in constructing top down 
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string duals. Most often, these phenomenological modifications are ad hoc constructions 
introduced to break the conformal symmetry and allow the introduction of the QCD 
scale. This was first investigated by Polchinski and Strassler [33], who truncated the 
AdS^ metric at some finite radius rmm- This approach is called the hard- wall model, 
and was further investigated by Erlich et al. in [M] . Rather than chop off a piece of the 
AdS space, Karch et al. [35 showed that the introduction of a non-trivial dilaton field 
(j){r) oc 1 /r^ also breaks the conformal symmetry. This latter approach is called the soft 
wall model, and is an improvement because the resulting meson spectrum displays the 
desired phenomenological behavior. One can also examine warped AdS metrics, which 
often have phenomenologically interesting thermodynamics. Some notable examples 
are the works of Gubser |36| 13 7j . Andreev [38] and Kajantie [39]. Finally, the work of 
Gursoy et al. [lU] HH US] 01] attempts to incorporate many of the features of QCD 
(e.g. running coupling) and are probably the most sophisticated examples of models 
built using the bottom up approach. 

2.2.3 Universality Approach 

Finally, one may hope that many or all strongly coupled theories share some common 
features. If in the course of working with other strongly coupled theories one discovers 
universal behavior which is model independent, one could then conjecture that the 
behavior would also be applicable to strongly coupled QCD. 

This approach has been especially fruitful in the case of the shear viscosity, which 
is a transport coefficient necessary for the hydrodynamic description of the quark-gluon 
plasma (more details will be given in Chapter [3|. Indeed, Kovtun, Son and Starinets 
showed that for a wide variety of gravity dual theories, the ratio of the shear viscosity 
rj to the entropy density s has the universal value 1/4tt. It was later shown that this 
relation holds for all theories of Einstein gravity |45| H6] (assuming the dual gauge theory 
is infinitely strongly coupled). It is quite remarkable that this result holds for both 
conformal and non-conformal theories, and is independent of the number of dimensions. 

One outgrowth of this universal behavior is the famous KSS bound conjecture [57] 
which states that for all physical substances. 



r]/s > l/47r. 



(2.20) 
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This has observable consequences for the plasma created at RHIC, even though the 
gravitational dual to QCD is not currently known. Hydrodynamic simulations of the 
quark-gluon plasma which attempt to fit the data from RHIC indicate that the QGP 
nearly saturates this bound |48| H9l [50l [51]. It is thus desirable to find other such 
universal behavior from gauge/gravity duality in hopes that there may be implications 
for heavy ion collisions at RHIC and at the Large Hadron Collider (LHC). 

The central results of this thesis employ this universality approach. In Chapters |4] 
and [5] we present the results for second order hydrodynamics for a wide variety of gravity 
duals. While these theories are not especially similar to QCD, one can examine the 
results to see if any other universal behavior is manifest. Some of the results presented 
in the aforementioned chapters could be applied to certain phenomenological bottom 
up models. 

2.3 Applications of the duality 

The idea that a strongly coupled gauge theory has a dual description in terms of an extra- 
dimensional gravitational theory has many practical applications. We will now review 
some of these applications. The application to hydrodynamic descriptions of the quark 
gluon plasma is the central focus of this thesis, and thus it will be presented in more 
detail in the next chapter. We will not dwell unnecessarily on the details of the other 
applications considered in the remainder of this chapter. Instead, a basic conceptual 
understanding of the relevant approaches is given, along with relevant references. 

2.3.1 Thermodynamics 

One can examine a strongly coupled gauge theory at finite temperature using gauge /gravity 
duality. In order to do so, one has to introduce a temperature into the gravitational back- 
ground. The obvious way to do this is to include a black brane (an extra-dimensional 
generalization of a black hole) into the background. Black holes have well defined 
thermodynamic properties, such as entropy and temperature; these thermodynamic 
properties are dual to the thermodynamics of the gauge theory. 

One gravitational background which is often used is the usual AdS^ x with the 
addition of a black brane horizon. The metric (in the near horizon limit r <^ L) takes 



the form 
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„2 t2 

ds^ = -^(-f{r)dt'^ + dx'^) + -—^dr^ + L^dnl (2.21) 
f{r) = 1-^. (2.22) 



Notice that this reduces to (2.4) if /(r) = 1. The position ro is the horizon; the time 
component of the metric vanishes at the horizon, while the r component diverges. We 
will often refer to this metric as the Schwarzs child AdS black hole, or simply as the 
SAdS metric. 

Rather than restrict ourselves to this metric, we will begin to use a more general 
black brane metric as originally considered by [47], 

ds"^ = gtt{r)dt^ + gxxir)dxidx^ + grr{r)dr^ . (2.23) 

Here, r is the coordinate denoting the extra dimension, and we assume r e (rg, oo). The 
index i runs over all spatial coordinates; let us denote the number of spatial coordinates 
by p. So that i = 1,2. ..p. The metric components are assumed only to be functions of 
r. There is p dimensional isometry along the spatial directions. 

The position of a horizon at r = ro is assumed. We define this position as the 
place where gu vanishes. If there is more than one such position, we take rg to be the 
maximum of these. Furthermore, we will need to specify the behavior of the metric 
components near the horizon: 

9tt{r^ro) ^ -7o(r-ro) + O(r-ro)2, (2.24) 

grrir^ro) ^ -^ + 0{1), (2.25) 
r - ro 

gxx{r ^ ro) gxx{ro) + 0{r - ro). (2.26) 
The quantities 70, 7r and ^^^(ro) are independent of r, but can in principle depend on 

The subject of black hole thermodynamics was pioneered by Bekenstein and Hawking 
|52^ [53l . The Bekenstein-Hawking entropy of a black hole for any theory of Einstein 
gravity is given by 

5=^, (2.27) 
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where is the area of the horizon. This quantity can be calculated from 



Ah = j V^dV (2.28) 

where "fij is the induced metric on the horizon, and dV is the induced volume element. 
These can be straightforwardly found from our metric by setting r = tq and dr = 0. 
Then, for our metric, 

y^^=m^- (2.29) 
There are many ways to derive the Hawking temperature from the metric. The 
approach listed below is based on the presentation of Zee [55^. Let us examine our 
metric in the near horizon regime. We will only need the t and r components for this 
argument: 

ds%H = -7o(r - ro)dt'^ + ^^dr^ + ... (2.30) 

r - ro 

Now, define a coordinate p so that 



P 



Then the metric is 



2V7r(r-ro), (2.31) 

dp' = (2.32) 

(r - ro) 



ds%H = -^p'^dt' + dp'' + ... (2.33) 



We now switch to an imaginary time variable r = it: 

ds%H= ^P^d^'^ + dp^ + - (2.34) 

In these coordinates, this piece appears to be written in polar coordinates, where p is 
the radial variable, and r is an angular variable. Examine the arc length a of a circle 
at constant p = po, and denote the period of the angular variable by /?. The arc length 
is 

a' = ]^pIp' (2.35) 

The distance a should be equal to the circumference of the circle, unless there is a 
conical singularity at the origin. We desire regularity at the origin, so 

(2vrpo)' = ]^pIP', (2.36) 
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Finally, as is customary in the case of finite temperature field theory, one equates the 
period of Euclidean time with the inverse temperature 

p-^ = T, (2.38) 

which gives us an expression for the Hawking temperature 

T=^,[^. (2.39) 
47r Y 7r 

Once we know the temperature and the entropy, it is easy to determine other quan- 
tities using thermodynamic identities. Throughout this thesis, we focus only on systems 
with zero chemical potential. In this case, the relations 

s{T) = (2.40) 
Ts{T) = e{T) + P{T) (2.41) 

allow us to determine the pressure P and the energy density e. Another quantity of 
interest is the speed of sound Vs in a thermally equilibrated medium: 

vl . (2.42) 

Again, in the case of zero chemical potential, we have 

" ~dJ/dT " Tds/dT- ^^-^^^ 
Generically, the temperature is a function of tq, and thus we can instead write 

2 ^ s(ro)((ir(ro)/dro) 

' T{ro){ds{ro)/dro)- ^ ' ^ 



Substituting in the previous relations for the temperature and the entropy (2.29), (2.39), 
we have 

1 -I^Logl^,] 
vl = - ^'-^^"^ . (2.45) 

This is a simple way to determine the speed of sound from the metric. We shall see 
in the next chapter that one can also compute the speed of sound by determining the 
hydrodynamic dispersion relation. This latter approach, which is much more compli- 
cated than the formula above, is of greater utility since it allows one to determine other 
transport coefficients of interest as well. 
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These thermodynamic tools are indispensable when attempting to study a strongly 
coupled gauge theory using gauge/gravity duality. As mentioned previously, these for- 
mulas for the entropy, temperature and speed of sound are used to study the thermo- 
dynamic properties of phenomenological models in j36l [371 [381 EHl Hll US] 

In addition to the thermodynamic aspects above, methods exist for examining phase 
transitions within the context of gauge/gravity duality. Specifically, the gravity dual 
of a phase transition is the Hawking-Page transition which occurs if more than one 
gravitational background satisfies the equations of motion derived from the action. In 
this case, one must examine the conditions under which each background is energetically 
favorable over the other(s). If at a critical temperature, one background is suddenly 
favored over the existing background, a transition occurs. This transition can be equated 
to the phase transition in the dual field theory. This methodology has led to holographic 
predictions of the deconfinement temperature \56\ l57] . 

2.3.2 Hadronic physics 

Much work has been done in describing bound hadronic states within the context of 
gauge/gravity duality. In the top down approach, as mentioned above, the introduction 
of flavor branes is necessary in order in include fundamental matter. Strings which begin 
and end on the D7 branes can be thought of as a quark/anti-quark pair. (Since the 
strings have no endpoints on the N D3 branes, they have no color indices, and are thus 
color singlets). Thus, excitations of this system correspond to the different hadronic 
states. There is a large literature on mesons in AdS/CFT, and an excellent review has 
recently been written 



One can also describe hadronic physics using the bottom up approach. This is 
done by introducing appropriate fields into the existing AdS^ background. The bulk 
fields are necessarily dual to operators in the field theory. The normalizable solutions 
to the equations of motion for these fields correspond to the hadronic states. For 
example, consider the left-handed current operator in QCD, qLl^t"'QL, where ql is 
a quark operator, and 7^, and are the usual Dirac matrices and group generators 
respectively. A suitable five dimensional field dual to this operator must be a gauge field 
j4^^°^ due to the Lorentz structure of the operator |^ Similarly, a right-handed gauge 

^ The mass of the dual field is set by the rules of the AdS/CFT correspondence. Here, m — [9lll0j. 
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field which is dual to the right-handed quark current operator is introduced. These 
two gauge fields can be combined into vector and axial fields '•"^ oc ^""^ + A^^""^ 
and A'^^") oc A^^""^ — A'^^"'\ One then constructs a Lagrangian involving these fields, 
and solves the equations of motion for the gauge field excitations. The normalizable 
solutions for the excitations of the vector field '•"^ are interpreted as the vector mesons 
in QCD. Similarly, the solutions for the axial field A^^^""^ correspond to axial mesons. 
Other techniques exist which allow for the calculation of other observables such as form 
factors and decay constants. 

These bottom up models are ad hoc. The backgrounds are often imposed by hand 
and not dynamically generated as the solution to any equations of motion; any back 
reaction which the fields have on the metric is usually neglected. The backgrounds are 
also not embedded within a high energy theory such as string theory. Despite these 
shortcomings, the phenomenological accuracy of such models is often close to 10%. The 
interested reader is referred to the review already mentioned [58] , as well as the primary 

works [SI Eg SD ESI EQl EU [62]. 

The behavior of mesons at finite temperature has also been investigated using 
gauge/gravity duality. By examining different configurations of the D7 flavor branes, 
one finds an interesting phenomenon related to the disassociation of mesons in a plasma 
at finite temperature. The idea is that the flavor branes may or may not end on the 
black hole. To determine which setup is energetically favorable, one needs to examine 
the free energy of each configuration. It has been shown that there exists a critical 
temperature above which the flavor branes end on the horizon (black hole embedding), 
and below which the flavor branes do not end on the horizon {Minkowski embedding). 
The spectrum of the mesons is very different in each of these two cases. In the low 
temperature phase, the spectrum is discrete and the mesons are stable; in the high tem- 
perature phase the spectrum is continuous and there is no mass gap. Thus, one observes 
a phase transition for the fundamental matter, which gives a holographic prediction for 
the temperature at which mesons disassociate in the plasma. This phenomenon was 
noticed in [631 El] , which built upon the earlier works [651 ES] . 
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2.3.3 Jet quenching and energy loss 

A jet is a large collection of particles which arise from the hadronization of a freely 
moving quark. Usually two jets are produced from a quark/anti-quark pair, and these 
jets are both seen in the detector. At RHIC, often only one jet is seen. The physical 
interpretation is that the missing jet lost energy while traversing the strongly coupled 
medium in order to reach the detector on the other side. This phenomenon is not 
observed in proton- proton collisions because a large region of quark gluon plasma is not 
created. 

We can learn a lot about a medium by examining how energetic particles traverse it. 
For this reason, it is desirable to examine how energetic quarks lose energy when travers- 
ing a strongly coupled plasma. This subject also helps us examine the phenomenon of 
jet quenching in heavy ion collisions at RHIC. From the point of view of gauge/gravity 
duality, this problem is usually approached in top down models. A quark, as mentioned 
previously, is the endpoint of a string on one of the flavor D7 branes. One can study 
a single quark by having one endpoint on the flavor brane, and the other on the color 
D3 brane (at the horizon). Alternatively, one can study a quark/anti-quark pair by 
having both ends of the string on the D7 brane. In either of these cases, the string 
connecting the two endpoints is described classically by the Nambu-Goto action. One 
can then examine the dynamics of this classical string in various situations (e.g. freely 
moving in one spatial dimension, falling toward the horizon, or under the influence of a 
constant force) . The dynamics of the string and its corresponding energy tell us about 
the dynamics of the quark(s) in the dual field theory. 

The pioneering papers in this regard are those of |67| I68j , as well as |69] which uses 
another approach; this latter work allows the calculation of a transport coefficient g, 
called the jet-quenching parameter using Wilson loops. These original papers have been 
extended in numerous ways, and this is still an active area of research. We refer the 
reader to the previously mentioned review [16] for a more complete set of references. 

2.3.4 Out of equilibrium and real time dynamics 

So far, we have only discussed the thermally equilibrated medium (QGP) that is suppos- 
edly created in heavy-ion collision. Of course, this grossly simplifies the problem, and in 
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order to have a full description of the collision, other effects must be taken into account. 
For example, we have considered an infinitely extended thermal medium whereas in a 
heavy ion collision the region has finite extent and expands after it is created. The 
introduction of flow into the gravity picture originated with |70| [7T] . 

Before and immediately after the collision, the system is not in thermal equilibrium, 
and the method by which it quickly equilibrates is not presently understood. From 
the point of view of the gravity dual, thermalization can be examined by perturbing 
a black hole background, and examining how long it takes to return to equilibrium. 
This approach has led to estimates of the thermalization time for a heavy ion collision. 
Typical estimates of this time are in the range of 0.3 - 0.5 fm/c, (see for example [72]). 
In addition to the thermalization time, one may also be interested in the isotropization 
time scale, and it can be accessed using similar methods |73| . 

2.4 Summary 

In this chapter, we have attempted to present the basics of the foundations of gauge/gravity 
duality and AdS/CFT. We have also explained three main philosophical approaches 
(top down, bottom up and universality) which attempt to adapt the existing AdS/CFT 
correspondence to the real world theory of QCD. Finally, we presented some example 
applications (and relevant references) of the duality to QCD and quark-gluon plasma 
phenomenology. 

The literature on this subject is vast and grows every day. This chapter is not meant 
to be comprehensive, but can serve as a starting point for the interested reader. This 
chapter also serves to present a context for the central results of this thesis. We have 
not yet discussed the application of gauge/gravity duality to hydrodynamics; this topic 
will be the focus of the remaining chapters. 



Chapter 3 



Hydrodynamics from 
Gauge/Gravity Duality 

This chapter is an introduction to the methods used in the central calculations of this 
thesis. We first review the theory of viscous hydrodynamics, and thus explain the origin 
of the transport coefficients; one of the main goals of this thesis is to calculate such 
coefficients. 

Secondly, we explain the methodology for calculating hydrodynamic dispersion rela- 
tions within the context of gauge / gravity duality. These dispersion relations give access 
to the transport coefficients. The necessary linearized Einstein equations are presented 
in full generality here. These equations will be solved in specific cases in the remaining 
chapters. 

3.1 Hydrodynamics 

In this section, we review the theory of hydrodynamics. To this end, we discuss how to 
construct the effective energy momentum tensor and the relevant normal modes which 
result from perturbations of this tensor. The discussion here follows excellent reviews 
presented in [75l ESI [771 ESI |I9] . 
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3.1.1 Hydrodynamic energy momentum tensor 

Hydrodynamics is best thought of as an effective theory which describes a thermal 
fluid when the length and time scales of interest are much longer than any relevant 
microscopic scale. In that case, one is able to smooth over the microscopic physics and 
instead use a perfect fluid description with viscous corrections. To lowest order, the 
energy momentum tensor takes the form of a perfect fluid 

r[^^ = -Pr/'^'^ + ^,V(e + P). (3.1) 

The subscript "PF" stands for "perfect fluid". Here, P is the pressure, e is the energy 
density, and is the fluid 4 velocity 

u'' = {jc,jv), (3.2) 

which takes the form (1,0, 0, 0) in the fluid rest frame. The flat space Minkowski metric 
is denoted by r/^i,; the fluids we consider will always be in flat space, while the extra 
dimensional gravity dual will have some non-trivial curvature. In general our conven- 
tions are those of Weinberg [76]. Our metric signature is ijf^^ = diag( 1,-1, -1,-1) in four 
dimensions, with an obvious generalization in higher dimensions. Greek indices indicate 
both space and time coordinates, while Latin indices (i, j...) indicate spatial coordinates 
only. 

One can proceed to add extra terms proportional to derivatives of the fluid velocity. 
In order to be in the hydrodynamic regime, the length and time scales of interest 
must be much larger than the microscopic scale, which is the inverse temperature T~^. 
Because the fluid variations are assumed to happen on long length and time scales, 
any these derivative terms are small corrections to the perfect fluid case. First order 
hydrodynamics comes from adding terms which contain at most one derivative. Let us 
now write these corrections in the local fluid rest frame where n* = 0, though in general 
derivatives of these quantities may not vanish. We define such a rest frame as the one 
where there is no flow of energy or momentum in the fluid; this is the so called Landau 
formulation. In such a frame, the fluid is at rest, so T"^ is still the energy density, 
and T*^* must vanish because there is no flow of momentum energy and momentum. 
Thus, the derivative terms can only show up in the spatial components of the energy 
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momentum tensor: 



(3.3) 



The viscous correction AT*-' is symmetric in i and j, and must be composed out of 
derivatives of the fluid velocity. Clearly, the two relevant terms which one can write 
are d'^u^ + d^u^ and rf'^d}^u^. Thus, the most general form of the first order derivative 
corrections is 



This equation is sometimes referred to as the constitutive relation. Here, we have intro- 
duced the variable p which denotes the number of spatial dimensions. We have written 
T*-' as above so that the first term which is proportional to r] vanishes under the trace 
operation. The coefficients r/ (shear viscosity) and Q (bulk viscosity) are the low en- 
ergy constants of the effective theory, and are referred to as transport coefficients. One 
should take care not to confuse the shear viscosity with the Minkowski metric. 

Physically, these coefficients describe a fluid's resistance to flow under stress. Shear 
viscosity is relevant for applied shear stress. As an example of shear stress, consider the 
following setup: lower a small, solid cylinder into a fluid and apply an external torque 
to the cylinder so that it rotates with a constant angular velocity. Naturally, the fluid 
very near to the cylinder will move with a greater velocity than the fluid further away. 
This velocity gradient will vary with the shear viscosity of the fluid. In fact, a similar 
technique is used to measure the shear viscosity of laboratory fluids, by measuring the 
amount of torque required to reach a certain angular velocity. 

Bulk viscosity is relevant for applied volume stress. To apply such stress, one should 
take a fluid element and compress it to a smaller volume (or expand it to a larger 
volume) without changing the shape of the fluid element. Bulk viscosity is thus only 
relevant for fluids which are compressible. 

Clearly, we are not in a position to directly measure these coefficients for the QGP, 
since it exists for an exremely short time. Instead, one must infer transport properties 
of the fluid from the multitude of hadrons which reach the detector after a heavy ion 
collision. Theoretically, these transport coefficients can in principle be calculated from 
the microscopic physics, but if the ffuid is strongly coupled, perturbative methods used 
to compute these quantities fail. Fortunately, gauge/gravity duality allows one to gain 




(3.4) 



25 

information about these coefficients. We will explain this how this is done in Sec. 3.2 

If the fluid under consideration has a conserved charge, one can also construct a 
similar relation for the current as an expansion in derivatives. Such an expansion will 
necessarily come with its own transport coefficients (e.g. the thermal conductivity k). 
Throughout this thesis we avoid any such complications; all theories under consideration 
here will have no chemical potential, and thus no conserved charge. In this case, the 



energy momentum tensor (3.4) is all that is necessary. 



3.1.2 Hydrodynamic modes 

Let us consider a fluid's response to hydrodynamic perturbations. To this end, we 
introduce fluctuations of the energy momentum tensor 6T^'^ and examine the equations 
of motion for such perturbations. We assume that all perturbations are of the plane 
wave type, e*^'^^""'*^ where q is the momentum and w is the energy of the perturbation. 
Without loss of generality, we have chosen our coordinate system so that z = Xp points in 
the direction of the momentum. By assumption, we require the energy and momentum 
of these perturbations to be small compared with the temperature w,q <^T (otherwise 
we are no longer in the hydrodynamic regime). 

Expanding T^'^ to linear order in these perturbations gives 

5r°^ = _5P^0m + (5e + 5P)n%^' + (e + P) (5n%^ + ^i°5u^) (3.5) 
= -dPri^f" + {5e + 5P)uf' + (e + P) {6u%'' + Su'') (3.6) 



and 



= -virtue + rj (^d'Su^ + d^5u' - '^r]'^dk6u^^ + Crf^ dk5u^ . (3.7) 



Here, we have introduced the speed of sound Vg defined as 

vl^^-f. (3.8) 
oe 

The equations of motion for such fluctuations can be found from the conservation 
of the energy momentum tensor 



5^ {ST^") = 0. 



(3.9) 
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First, let us consider the case of = 1, recalling the assumed plane wave style depen- 
dence of the perturbations. We have 

dtST^^ + d,5T^P = (3.10) 



Substituting results from (3.6) and (3.7), 

5T^^ = ^e + P)5u^ (3.12) 

8T^P = 7j {d'5{u^)) = -{iq)7j6u^, (3.13) 

we have 

{-iw){e + P)6u^ + q'^r]6u^ =0, (3.14) 
which requires a dispersion relation of the form 

w{qUe.r = -i^^o" + 0{q^). (3.15) 
e + P 

The corrections of order q^ to the dispersion relation would come from higher order 
dissipative corrections. Since we have only written the energy momentum tensor to first 
order in derivatives, our dispersion relation is only valid to the leading order in q. This 
is the shear mode perturbation; its dispersion relation depends on the shear viscosity. 
We derived it by assuming z/ = 1 in the conservation equation, but we would get the 
same relation if we used u = 2,3...p — 1. 



Next, consider the conservation equation (3.9) with u = and u = p. We then get 
two coupled equations 

- iw6T°° + iqST^P = 0, (3.16) 
-iwdT'^P + iqSTPP = 0. (3.17) 



One can now go back to the constitutive relations for JT^"^ (3.4) and substitute them 
here. Before doing, so, however, it is useful to note that when introducing these per- 
turbations, we desire preservation of the normalization of the velocity vector. This 
means 

S {rif^^u^u") = (3.18) 
ri^f^{2u'^5u'^) = (no summation) (3.19) 
6u° = (fluid rest frame). (3.20) 
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In the last line, we have specified to the local rest from where only is non-zero. Using 



this fact along with (3.6), (3.7), we can now write 



6TPP 



<5e, 

(e + P) SuP, 
Vg6e — iqSu^ 



C + 2 



p — 1 



Plugging these back into (3.17) gives 

-w (e + P)5vF + q{ v^Je - iqdu^ 
which combine to give 

5uP (e + P) 



w5e + q{e + P) Su^ 



2 I 2 2 

-w + q v„ 



iq^w 



e + P 

After solving for w{q) and expanding in powers of q we find 

.2 r 



0. 



w{q)sonnd = ±Vsq 



iq 



C + 2 



P 



7] 



+ 0{q'). 



(3.21) 
(3.22) 
(3.23) 



(3.24) 
(3.25) 



(3.26) 



(3.27) 



2{e + P) 

Our hydrodynamic dispersion relation is only accurate to first order in the dissipative 
terms. In order to go to higher order in q, we would need to add terms with more than 
one derivative to the effective energy momentum tensor. 

This latter dispersion relation is the sound mode (or compressional mode). It corre- 
sponds to a perturbation moving with speed Vg relative to the fluid; the dissipation of 
the wave is controlled by a combination of the shear and bulk viscosities. 



3.1.3 Second order (causal) hydrodynamics 

In the preceding section, we discussed first order hydrodynamics; in this section we 
discuss the next hydrodynamic order. This subject was first approached by Miiller |80j 
and later Israel and Stewart |8H |82] . It is desirable to extend the theory of hydrody- 
namics to the next order, because the first order theory has problems with causality 
|83j . Formally, the issues with causality have been shown to exist only for modes which 
are outside the hydrodynamic regime |84| [85| \8E\ [57] . but from a practical standpoint. 
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such issues are unacceptable when attempting to do numerical simulations |86| [871 EH] • 
These are the issues which the Israel-Stewart formulation attempts to resolve, and for 
this reason the approach is sometimes called causal hydrodynamics. 

The examination of the next hydrodynamic order is also attractive from a theoretical 
standpoint. Gauge/gravity duality has led to the important observation that in all 
Einstein gravity duals, the ratio of the shear viscosity to entropy density r]/s takes 
on the universal value I/Att (this will be discussed in more detail in the subsequent 
chapters). Given this success, it is interesting to inquire whether other such universal 
relations exist at the next hydrodynamic order. 

Israel introduced five new transport coefficients that appear in the hydrodynamic 
expansion of the energy momentum tensor. In what follows, we use the same notations 
and conventions as [SH] • Three of these five transport coefficients are relaxation times 
associated with the diffusive, shear, and sound mode, and are denoted by (tj, T,r, ''"n) 
respectively. There are two other transport coefficients which are related to coupling 
between the different modes, ao,cei. 

We will not write down the full expression for the second order energy momentum 
tensor in Israel-Stewart theory, but refer the interested reader to the original works 
mentioned above. For our purposes, we need the expression for the shear and sound 
mode dispersion relations in terms of the transport coefficients. These relations were 
worked out by Natsuume et al. [89^ for the case of a background without a conserved 
charge. (As mentioned previously, this is the case which we focus on here. A more 
complete list of assumptions regarding this dispersion relation can be found in |89j ) . 
The results for the sound mode are 



w{q)sound 



wiq + W2q^ + w^q^ + 0{q'^) 



(3.28) 




(3.29) 



29 

with wi and W2 unchanged from ( 3.27[ )p] Furthermore, 

w{(l)she.r = -ij^q^ - i (^-^^ (r^ + A)(?4 + 0(g6). (3.30) 

A few words are in order about the parameter which we have denoted as A above. 
Suppose that we are after the next order correction to the shear mode dispersion relation. 



Let us return to the conservation equation for the shear mode (3.13) 



- u;<5r°i + q6T^P = 0. (3.31) 

In the rest frame which we consider, 6T^^ ~ C'(l), it does not depend on q. If we 
are working within the context of second order hydrodynamics, 5T^p can have terms 
proportional to one or two derivatives. In general, it will have the form 

6T^P = Aq + Bw + Cw^ + Dwq + Eq^, (3.32) 

where A, B, C... are functions of the fluid four velocity, but are independent of q. The 
conservation equation becomes 

- W(5r°i + + Bwq + Cw^q + Dwq^ + Eq^ = 0. (3.33) 

Expanding w{q) in even powers of g, one finds that the term which is relevant for 
determining the next correction to the dispersion relation is the term containing D, 
since this is the only piece which has a term proportional to q^. However, it is clear 
that there is another way to produce a term which is proportional to g^, and that is 
if T^P ~ g3_ Clearly, such a term could occur from three spatial derivatives, and is 
thus only possible within the context of third order hydrodynamics. Hence, we have 
shown that the correction to the shear mode dispersion relation is dependent on the 
third order hydrodynamic expansion of the energy momentum tensor. This fact was 
first realized in |90j. One might be surprised by the fact that in order to consistently 
compute the subleading term in the hydrodynamic dispersion relation, one needs to use 
third order hydrodynamics. This is simply a consequence of the fact that the to lowest 
order 'u;(g)shear oc g^. 



^ Note that there are actually two solutions here, corresponding to the plus/minus sign in the 
formulas for wi and U13. These two solutions only differ by the relative direction of the momentum q. 
Throughout this thesis, we will always assume that our coordinate system is chosen so that Vs > 0. 
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Currently, there is no formulation for third order hydrodynamics. We have intro- 



duced the parameter A which appears in the shear mode dispersion relation (3.30) to 
denote the sum of all such contributions from third order hydrodynamics. 

3.2 Dispersion relations from gravity 

Observations of elliptic flow of the plasma created at RHIC indicate that the system 
exhibits collective motion characteristic of a thermally equilibrated system, and is thus 
describable by hydrodynamics [H |5l |6l [71 [911 [92l [93]. The transport coefficients men- 
tioned in the previous section are necessary input for hydrodynamical simulations of this 
matter. As outlined in Chapter |2] gauge/gravity duality allows one to map a strongly 
coupled gauge theory to a classical gravity theory in an extra dimension. Thus, there 
exist methods for determining transport coefficients, and the shear and sound mode 
dispersion relations within this framework. 

There are many ways to calculate the transport coefficients and dispersion relations 
of a strongly coupled gauge theory using an extra-dimensional gravity dual. One can 
compute correlation functions of the stress-energy tensor and use Kubo formulas or 
examine the poles of such correlators |94| [95l l96l [971 IHHl |99]. Alternatively, one can 
examine the behavior of the gravitational background under perturbations and deter- 
mine the dispersion relation for such perturbations by applying appropriate boundary 



conditions. Comparison with the expected dispersion relations from (3.151,(3.27) yields 
formulas for the transport coefficients [21 11001 HOI] . In addition, the black hole mem- 
brane paradigm has been employed to calculate the hydrodynamic properties of the 
stretched horizon of a black holej^ In many cases, the transport coefficients cal- 
culated on the stretched horizon coincide with the transport coefficients in the dual 
gauge theory [H [13 [HI [Ml [Ml [Mjjf] Recently, the work of [106j provides yet 

^ This is the idea that a black hole's influence on the outside world can be encoded in an effective 
membrane which lies just outside the horizon. This idea is similar to replacing a spherical mass dis- 
tribution with a point mass; provided one remains outside the event horizon, the membrane appears 
physically equivalent to the actual black hole. This effective membrane is sometimes called the stretched 
horizon and can be endowed with thermodynamic and hydrodynamic properties such as electrical con- 
ductivity and viscosity |102| . 

^ For an example of a situation where the stretched horizon transport coeflicients differ from those 
in the dual field theory, one can consider the bulk viscosity. The bulk viscosity on the stretched horizon 
is negative, whereas the bulk viscosity of the dual field theory is non-negative 29 • 
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another way to compute hydrodynamic transport coefficients (sometimes referred to as 
fluid/gravity correspondence), by deriving the equations of fluid dynamics directly from 
gravity. This work has proved quite influential, and has led to much subsequent research 

[Ml [Ml [Ml [nni Ell ttn] • 

In the past few years, much work has been done to extend previous analyses to 
second order hydrodynamics [H |89l ISOl IMl [Ml IMl [Ml ttHl [112]. Most of the 
work on second order hydrodynamics has focused on conformal theories. It is notable 
that a universal relation between second order hydrodynamic transport coefficients of a 
conformal theory was presented in |113] . though it is not known whether this relation 
still holds for non-conformal theories. 

In this thesis, we use the gravitational perturbation approach similar to |100| llOlj 
to compute the shear and sound mode dispersion relations for a special class of gravity 
duals. In this section we explain this methodology. We will always work with a black 
brane type metric (2.23) as discussed in Chapter [2j 



The idea behind this method is conceptually simple. Fluctuations in the strongly 
coupled plasma should be dual to fluctuations of the gravitational background. The 
latter are introduced by adding perturbations in the metric 

9fiu gfiu + hf,^ (3.34) 

and any matter fields present. For example, if a scalar field is present, one must introduce 
the fluctuation 

(/> ^ (/> + (5(/). (3.35) 

In later chapters, we will specify the matter supporting the metric as scalar in nature. 
For now, though, we will work with a general energy momentum tensor. It is important 
to note that the energy momentum tensor discussed here describes the matter which 
supports the gravity dual. It is not the same as the hydrodynamic energy momentum 
tensor discussed in the previous section. This latter quantity is defined on the (p+1) 
dimensional boundary of the bulk space-time. In short, one should take care to distin- 
guish the energy momentum tensor which supports the extra-dimensional metric, and 
the energy momentum tensor of the dual field theory. 

Returning now to our perturbations of the gravity background, one must take the 
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background Einstein equations 

Gg] = -87rGp+2rW, (3.36) 

and expand them to linear order in these perturbations. Symbohcally, we write 

G(,i) = -SttGp+sTW. (3.37) 

Here, G^y is the Einstein tensor and Gp+2 is the p + 2 dimensional Newton's constant. 
The superscripts (0) and (1) denote the order of the perturbation. Explicit expressions 
for the background and perturbed Einstein tensor are given in Appendix [A} (In addition, 
the background equations for the matter fields must also be expanded to first order in 
the perturbation). 

The resulting set of linearized equations can then be solved for the perturbations. 
The perturbations of the hydrodynamic energy momentum tensor were considered to 
be plane waves in the previous section. Hence, we should make the same assumption 
here, but there is an extra complication owing to the fact that the metric lives in one 
extra dimension. Thus, we make the ansatz 

V = VWe^^''""'*^ (3-38) 

and similarly for the matter perturbations. One must then solve the linearized Einstein 
equations for the functions /i^i^ perturbatively in q (since we are still working in the 
hydrodynamic regime w,q <^ T). Applying appropriate boundary conditions to these 
solutions will result in a dispersion relation w{q). One can then read off relations for the 
transport coefficients by comparing the resulting dispersion relations to those expected 
from the hydrodynamic fluid (3.15 ), p. 27 ). This allows one to write expressions for 



the quantities in the dispersion relation {wi,W2^ •••) in terms of the metric components 
idtty 9rr, dxx) and any background fields. Let us now discuss the individual steps of this 
method in more detail. 

3.2.1 Classification of hydrodynamic modes 



In Sec. 3.1.2, we described how perturbations of the hydrodynamic energy-momentum 
tensor can be decomposed into two normal modes: shear and sound. Since we assume 
that our system admits a dual gravitational description, we need to do the same thing 
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for the gravitational perturbations. Which components of the metric fluctuation h^i, 
correspond to shear /sound perturbations in the dual field theory? 

These components are determined by classifying the perturbations under the rota- 
tion group SO{p — 1) [13 [96 l 1101] . We have singled out the z direction as the direction 
of our momentum, but there is still symmetry in the remaining p— 1 spatial dimensions, 
and thus such a classification is possible. A general transformation matrix for a spatial 
rotation about the z-axis is given by 

/ 1 ... ^ 





U{9) 









Ai,iW 



Xp-i,i{e) 










Ap-i,p-i(e) 
1 

1 



(3.39) 



Here ^ denotes the row, and i' denotes the column. The indices are ordered such 
that fj, = 0,1, ...p— l,p,p + 1 corresponds to t, xi, X2..., Xp^i, z,r. The components 
A(^) describe how the various spatial components are rotated into one another. (As a 
simple illustration, in the case of three spatial dimensions, Ai^i = A2,2 = cos(0), and 
Ai^2 = — A2,i = sin(^)). Let us now examine how the components of /i^^ transform 
under this transformation 



For example, consider 



af3 



11°' h 



(3.40) 



(3.41) 



In what follows, we use the Latin indices (a, 6, c.) to denote the spatial coordinates 
xi...Xp-i. Then, it is clear that 



htr 


— htr 


(3.42) 


hzr 


— hzr 


(3.43) 




= hrr 


(3.44) 


har 




(3.45) 
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Notice that the components htr, hrr, hzr transform into themselves, and are thus referred 
to as scalars under this transformation. The components hbr transform with one factor 
of U, which is how a vector transforms. Continuing with other components, one finds 
that hzz,htz,htt belong to the scalar mode, while hta, and hza belong to the vector 
sector. The only unaccounted for elements are the hab components. 

The components of hat can be split into a trace and a trace- free part. How does 
the trace piece transform under the spatial rotation? Clearly, our rotation is a unitary 
transformation. Recalling that the trace of a matrix is unchanged under a unitary 
transformation, it is required that 

htt + '^haa + hzz + hrr = htt + ^ haa + hzz + Kr- (3.46) 
a a 

Since we've already shown that htt,hzz, and hrr are scalars under the transformation, 
we conclude that the trace of h^t is also a scalar. 

The remaining (trace-free) part of hab transforms as a rank two tensor. These com- 
ponents constitute a separate mode (the tensor mode), which we will not be concerned 
with in this thesis because the the resulting dispersion relation does not have a well 
defined hydrodynamic limit (see for example |89]). In full, we have the following tensor 
decomposition for the metric perturbations 

htt , Zla ^a.a , hzz , Kr , hz , hr , Kz Scalar mode (sound mode) 

hra, hta, hza vector mode (shear mode) (3.47) 

hab - ^ab:^ Ylc hcc tensor mode . 

3.2.2 Shear mode equations 

In the previous subsection, we found that to describe the shear mode, we should make 
only the following components of the metric perturbation non-zero: hta, hza, and hra- 



We now make the ansatz mentioned in the Sec. 3.2 and define the functions M{r) and 
N{r) as 

hta = /iia(r)e^(^^-"'*)^fcM(r)e^(«^-"'*), (3.48) 
hza = hza{r)e'^'''-'"'^ ^g..N{r)e'('''~^'l (3.49) 

There is always relativistic gauge freedom when one deals with gravitational perturba- 
tions. We are free to choose our coordinate system so that /i^r vanishes. This is referred 
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to as the radial gauge. It is important to note that this choice only partially fixes the 
gauge; later we will define a set of gauge invariant variables to deal with the remaining 
gauge freedom. 

In Appendix |Xj the relevant formulas are given which allow one to compute the 
Einstein equations to linear order in the metric perturbation. Using these formulas 
with only the shear perturbations non-vanishing results in three non-trivial Einstein 
equations]^ These equations come from the {ta), (za) and (ra) components of the 



linearized Einstein equations (A. 35). They can be written: 
1 



-9 



dr [^g"'g''g.xM'] - qg'' {qM + wN) = -167rGfc T/f) + r/*)M , (3.50) 



--dr [V^g^'^'N'] - wg'' {qM + wN) = -IGnGg, 



J(i) + J(o)^^ 



qg^'^'N' - wg^^M' 



(3.51) 
(3.52) 



Here, the prime denotes derivatives with respect to r. When not indicated otherwise, G 
is short for Of course, components of the background energy momentum tensor 



T^^^ can be written in terms of the metric components using the relations (A. 21 - A. 23 ) 



3.2.3 Sound mode equations 

In contrast to the shear mode, the relevant perturbations which must be considered for 
the sound mode are hu, htz, hzz, and haa- Again we have partially fixed the gauge by 
setting hnr = 0. We proceed as before by defining 



p-i 



1 ^ 

a=l 



hu{r) 

haa{r) 

hzz{r) 
htz{r) 



gtt{r)A{r), 
gUr)B{r), 

gxx{r)C{r), 
gtt{r)D{r), . 



(3.53) 

(3.54) 

(3.55) 
(3.56) 



* One may wonder why there are only three equations, when at first sight there are 3(p— 1) degrees 
of freedom in this channel. The equations for a = xi, for example, are identical to those with a = X2 
up to a change of variable. This is why we have defined M and A'^ as above, and not taken care to 
distinguish M\,M2, etc. 
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Ostensibly, there are 5+p Einstein equations. These come from the {tt), {zz), (rr), (tz), 



{rz), (tz) components of the Hnearized Einstein equation (A. 35), plus the {p — 1) equa- 
tions for the (aa) components of this equation. However, each of these latter equations 
is identical up to a change of variable. For example, label two of the a coordinates as 



X and y. Then the equation involving the {xx) component of (A. 35) is identical to the 



equation equation involving the (yy) component up to a replacement h^x — > hyy. We 
can add all of these p — 1 equations together and get a single equation which involves 
only the variable B defined above. Thus, there are seven Einstein equations in total. 
To simplify the presentation of these equations, we make the following definitions 



f{r) s -S„s"- 



X r) 



(3.57) 
(3.58) 



Then, the relevant seven Einstein equations (in the radial gauge) can be written 



fdr 



-99 

V7 



-{{p-l)B' + C') 



{p - l)g'5rr5""i? (3.59) 



dr 



yZ^^-^'l + ^dr[V^gg^^{{p-2)B' + C')] 



(3.60) 



9rrg''''\ 



911) 

A + {p-2)B + —D 



w 
J 



[ip-2)B + C] 



IQ^Ggrrgxx E + Tl^^B^ 



dr 



y^fg^'^A'] + -^dr[V^gg"'{p-l)B'] 



(3.61) 



{p - l)w^grrg''B = IQnGgrrgxx (T^r) + ^^of I , 



dr [\/=^/7^'] + (P - l)qwg''grrB = -IdnGg^^grr [Tf^) + Tfop) , (3.62) 



VL[{gxxr]A' + VL[gu{gxxr-']{{p-l)B' + G') 

- 2grr\w'g''[{p-l)B + G]+q^g^^ 



(3.63) 



2w 

A + {p-l)B + — D 



327rG(5rr)'T(7), 
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9 « 



-^((p-l)i? + C) 



q!D' = -I (l67rG5,,ff«r*[)) , (3.64) 



L V . ^^j + y^r [/^] + g(p - 1)5' = i (l6TTGgrrgxxT(,]j . (3.65) 



7^ 

These equations are the Einstein equations which involve the {tt), (aa), (zz), {tz 



(rr), {tr) and {rz) components of the (A. 35) respectively. We emphasize that these 



equations are valid for any matter distribution which is minimally coupled to the metric 
(in other words, they are valid for Einstein gravity, and any matter distribution). When 
a particular form of matter is chosen, the components of the energy momentum tensor 
on the right side of the above equations will have to be explicitly evaluated. 

Whatever matter fields are present will obey separate background equations. In- 
troducing perturbations into these equations will result in equations for the matter 
fluctuations (e.g. 5(j) for a scalar field, 5A^ for a gauge field, etc.). We will determine 
the form of these equations for scalar matter in Chapters |4] and [5} 

3.2.4 Gauge invariance 

When examining perturbations of a general relativistic background, it becomes difficult 
to tell the difference between physical perturbations, and those which could be inter- 
preted as a coordinate transformation. The ability to change ones coordinates at will, 
is referred to as relativistic gauge freedom. In particular, metric perturbations which 
are related to each other under the transformation 

V ^ V - Vj'^e. - y^uH^. (3.66) 

describe the same physics. Here is any vector, and V^''^ denotes the covariant deriva- 
tive with respect to the background metric (see Appendix [A|) . In order to preserve the 
assumed space-time dependence of our perturbations, the vector should also contain 
the same dependence 

i^. = Ury^'"~""^- (3-67) 

In order to address the issue of gauge freedom, one can take two approaches: either 
fix the gauge completely, or do the analysis in terms of gauge invariant variables. In 
this thesis, we use the latter approach. Hence, before we determine our dispersion 
relation, we will need to create gauge invariant combinations of the perturbations which 
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do not transform under (3.66). The particular choice of gauge invariant variables will 
be given explicitly in the later chapters, as they depend on the type of matter fields that 
are present. Gauge invariant variables have been used extensively within the context of 
cosmology (cf. |114] ): in the context of gauge/gravity duality, they were first introduced 
in |100] . In this thesis, we use Z to denote gauge invariant variables; sometimes we will 
also add an additional subscript (e.g. Zq, when it is necessary to make a distinction 
between several such variables. 

When constructing such gauge invariant combinations, one has two possible choices. 
The first option is to first solve for the perturbations h^u{r), and then construct gauge 
invariant combinations from these solutions. The second option is to first take combi- 
nations of the linearized Einstein equations and reduce them to a set of equations which 
depend only on the gauge invariant variables, and then solve these equations. In our 
analysis of the shear and sound modes, we will employ the latter approach. 

3.2.5 Boundary conditions 

After deriving the relevant equations for the perturbations, one must solve these equa- 
tions perturbatively (order by order) in the variable q by first expanding the pertur- 
bations in powers of q, and also inserting an ansatz for the dispersion relation w{q). 
We are ultimately interested in the dispersion relation (not the solutions for /i^;/(r) or 
Z{r)). One must apply appropriate boundary conditions on the solutions in order to 
determine w[q). 

Two obvious places to apply boundary conditions are at the horizon and at the 
boundary r — > oo. The correct boundary conditions that need to be applied on the 
gauge invariant variables were worked out in [100 . In this work, the authors constructed 
gauge invariant equations for the AdS^ background. By examining the behavior of the 
equations in the near horizon regime, one finds that the gauge invariant variable Z must 
behave as 

Z(r^ro) ~ (r - ro)" + 0(r - ro)°+^ (3.68) 

a = ±-^. (3.69) 
47rr ^ ' 

One interprets the positive/negative exponent by stating that the perturbation corre- 
sponds to either an outgoing/incoming wave. Because classically nothing can be emitted 
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from the horizon of a black brane, we should choose the mmus sign. This is referred to 
as the incoming wave boundary condition. It can be implemented by making the ansatz 



Z(r) = /(r)--/4-^y(r), 



(3.70) 



where the function Y(r) is regular at the horizon - it either vanishes there or approaches a 
constant. Another way to implement this condition is to take the logarithmic derivative 



of the above equation (using the notation defined in (3.58)) 



VL[Z{r) 



iw 
1^ 



VL[f{r)]+VL[Y{r)]. 



(3.71) 



Let us examine the leading behavior of this equation near the horizon, using the near 



horizon behavior of the black brane (2.24 - 2.26) and the fact that Y{r) is regular at 
the horizon. This leads to 



Z'{r ro){r - tq) 



Z{r ro). 



(3.72) 



A condition at the boundary (r oo) is also needed. The authors of |1UU] show 
that the correct boundary condition to apply is 



Z{r 



oo 



0. 



(3.73) 



The authors of the above mentioned paper argue this is the correct boundary condition 
because if one applies this condition, one correctly reproduces the dispersion relation 
that appears as a pole in correlation functions of the energy-momentum tensor. In other 
words, if one applies this boundary condition, one finds agreement between this method 
and other methods of computing the dispersion relation. 



3.2.6 Summary 

To summarize, in order to compute hydrodynamic dispersion relations using gauge/gravity 
duality, one must follow the following prescription. 



1. Depending on which hydrodynamic mode is being studied, one must introduce 
the relevant metric perturbations (3.47), as well as perturbations of any and all 
matter fields present. 
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With the relevant perturbations turned on, one must hnearize the background 
Einstein equations, and equations for the matter fields. The linearized Einstein 



equations for our particular black brane metric are given explicitly in (3.50 - 3.52 ) 



for the shear mode, and (3.59 - 3.65) for the sound mode. 

3. Either: 

• Reduce this set of equations down to equations which depend only on gauge 
invariant variables and then solve for the gauge invariant variables order by 
order in the momentum of the perturbation, q. 

Or: 

• Solve for the perturbations h^u{r) order by order in q, and construct gauge 
invariant combinations from these solutions. 

4. For the resulting solutions, apply the incoming wave boundary condition in the 



form of (3.70) or (3.72). Also apply a Dirichlet boundary condition at the bound- 
ary (r oo). 



In the next two chapters, we will illustrate this procedure for the shear and sound modes 
assuming the matter fields present are scalar fields. 



Chapter 4 

Shear mode 



4.1 Introduction 

In this chapter, we will apply the prescription outlined in the previous chapter to the 
shear mode. We assume that the matter which supports the black brane metric is scalar 
in nature. To this end, we must first explicitly compute the components of the energy 
momentum tensor to first order in the perturbations. Once this is done, we proceed to 
construct gauge invariant equations, solve them, and apply the appropriate boundary 
conditions to determine the shear mode dispersion relation. Shear mode perturbations 
are strongly over-damped, and the shear mode dispersion relation is an expansion in even 
powers of the momentum q. Let us parametrize the dispersion relation by introducing 
two real constants -D^ and Tshear- 

W{q)sheiii = -iDr^q^ (l + rehear -D.^ 9^ + •••) • (4.1) 

In what follows, we will determine and Tshcar in terms of the metric components 
gtt,gxx, and grr- Through comparison with the hydrodynamic expectations for the 



shear mode (3.30), the computation of allows one to derive a formula for the shear 
viscosity t], and the result for Tghear provides information the relating second order 
transport coefficient to the correction from third order hydrodynamics A which was 
introduced in section 13.1.31 

The formula presented here for Z)^ was first derived by Kovtun, Son and Starinets in 
j47j using the membrane paradigm, and later in [115j using techniques similar to those 
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we employ here. The formula for Tghcar was first derived by Kapusta and Springer in 
m again using the membrane paradigm. This chapter is based on this previous work, 
though here we do not make any reference to the membrane paradigm. As far as we are 
aware, this is the first instance in the literature where the formula for Tghear is derived 
without any reference to the black hole membrane paradigm. 

Finally, once we have the formulas for L)^ and Tshcar in hand, we will consider a few 
applications of these formulae to specific gravitational backgrounds. 



4.2 Scalar matter 

4.2.1 Background equations 

Let us assume that the matter supporting the metric is a set of n minimally coupled 
scalar fields. In other words, we assume the action is of the form 



1 



levrG 



p+2 



dP+^xV^ {R - ^d^^Pkd^cPk - U{(t>i,(t>2:4n)) ■ (4.2) 



We assume a summation over the repeated index k, and we use superscripts (0) and (1) 
to denote background quantities and quantities that are first order in the perturbation 
respectively. Then, the background equations are 

Gji) = -SttGp+sTW (4.3) 
O^'Uk = (4.4) 



Here O^^^ is defined as in the appendix (cf. A.18). We ignore any subtleties regarding 



boundary terms that come from integration by parts. Such terms can be taken care of 
by adding additional boundary terms to the action. 

The energy-momentum tensor derived from the action is 

87rGp+2T^? = l{d^(l>kdu(l>k-9t.uC^), (4.5) 
= ^dx^kd^<pk + Ui(l)i,(l)2.:M. (4.6) 

Because the background metric only depends on the extra dimensional coordinate r, if 
n = 1 it is clear that the single field (j) must also only a function of r. In the case of 
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multiple scalar fields (n > 1), it might be possible to have fields which depend on the 
other coordinates, provided that all such dependence cancels out in the combination 
that appears in xj^^ . We will not consider such special cases, and will always assume 
that the scalar fields only depend on r. 

One can write the fields (pk (f) in terms of the metric components by noting that 

where the F functions are defined in the Appendix lAj Furthermore, 



n 



(4.8) 



k=l 

Thus we have the relation 

n 

J2 4>'k{r? = 2grr [Ft{r) - F,(r)] . (4.10) 

k=l 

It is also noteworthy that this background has the special property Ft = as can be 
seen by considering 

5«g(?-5^^gW = -8vrG,+2 (^"T;?^ - rtS)) (4.11) 
Ft{r)-F,{r) = £^(r) - £^(r) = 0. (4.12) 

Because of this fact, and the general theorem given in all backgrounds we 

consider saturate the conjectured shear viscosity bound: rj/s = 1/4tt. 

4.2.2 A note on phenomenological model building 

Let us make a small digression and examine implications of these background equa- 
tions for phenomenological model building. As mentioned in Chapter |2] the bottom up 
approach (or AdS/QCD) often employs metrics and scalar fields which are chosen to 
reproduce some essential features of QCD. Often, these metrics are imposed by hand, 
and are not dynamically generated from any matter distribution. 
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One such popular background is the soft-wall model of [35]. At zero temperature, 
this five dimensional model has a (string frame) metric and background scalar field 
(dilaton) of the form 

7-2 1,2 
'^■Sstring = j2 (~'^*^ + dxidx'^) + ^dr^ , (4-13) 

(/>(r) = ^. (4.14) 

Here, L is the usual AdS curvature radius and c is a dimensionful constant which allows 
for conformal symmetry breaking and the introduction of the QCD scale. In |116| . it 
was shown that such a background can be realized dynamically with the addition of a 
second scalar field x(^)) ^"^^ appropriate scalar potential V{(j),'x). 

These results have not yet been generalized to finite temperature. Here we will 
show that an obvious generalization of this metric to finite temperature is not possible. 
Because the string frame metric is exactly AdS^^ a naive expectation is that at finite 
temperature, the metric should take the form of the Schwarzschild AdS black hole. 
Thus, we might expect a finite temperature generalization of the metric to be 

^^•Sstring = j2 {-fir)dt'^ + dxidx') + ^^^^^^ (ir^ , (4.15) 

fir) = 1-^. (4.16) 

The dilaton is a scalar field which is coupled to gravity by a term e'"^R in the (string 
frame) action where 6 is a constant, and R is the Ricci scalar. Thus, in the string 



frame, the action is not of the form (4.2). However, it is easy to make a conformal 



transformation of the metric which brings the action into the form (4.2) (cf. Appendix 
G of |117j ). This is called the Einstein frame; in this frame, the metric takes the form 

yEinstem ^ ystrmg- K^-'-') 

a is a constant which depends on b and the number of dimensions of the theory; its 
value is unimportant for our purposes. 

With the Einstein frame metric in hand, we can use the results of the previous 
section. For example, if this metric is generated by scalar fields, it must obey the 



constraint (4.12). An explicit evaluation of this equation yields 



^2,3 - 0- (4-18) 
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Clearly, this is inconsistent as long as a is non-zero and the dilaton (j) has some non-trivial 



profile. The conclusion is that the metric (4.15) cannot be generated by scalar fields 



alone. One can phrase this another way by stating that one must modify the metric 
away from the SAdS form in order to generate a finite temperature soft wall model 
using scalar fields alone. The reason why there is no problem at zero temperature is 
that one can arrive at the zero temperature AdS^ metric by taking the limit tq — > 0. In 



this limit, the equation (4.18) is satisfied and there is no inconsistency. 



4.3 Shear mode equations for scalar matter 

Let us now move beyond the background equations and examine the relevant shear 
mode perturbations. In order to make use of the general Einstein equations given in 



the previous chapter (3.50 - 3.52 ), one must determine the components of the perturbed 



energy momentum tensor T^^^. For the background, we have 

87rGp_,2T,'^0) = I {9'"'9''^d^<t>kdp4>k - 9^"" ^) ■ (4.19) 



(0) - 2 

Including the perturbations g^^ g^^ — h^^ , and (pk 4>k + S4>k, to linear order we 
have, 

87rGp+2Tl{] = l{g''''g''f'[da{6^k)dp<t>k + da(t>kdp{6(Pk)] 

- h^'^g^^do.(t>kdp<i)k - g^''h^''da(t>kdp<Pk 

- g^''{8C^) + h^''C^]. (4.20) 

With the choice of the radial gauge h^^ = 0, the second line vanishes due to the fact that 
the background scalar fields only depend on r. For the shear mode, we need to evaluate 
T*", T^", and T'~°. For linear response, it is natural to assume that the perturbed scalar 
fields have the same space-time dependence as the metric perturbations 

50fc = <5(^fc(r)e^(^^-'"*). (4.21) 

Then, these three components of the energy momentum tensor are straightforwardly 
given as 

BTrGp+sT^f) = \h'^C^, (4.22) 



2 ^<f>1 



0. 



Noting that for the background, 



1 „xx n 



off '-'t'l 



we have the relations 



rpta 

^(1) 

rpza 

^(1) 



(0) 



-M{r)Tf' 



(0) 



(0) 

(0)- 
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(4.23) 
(4.24) 

(4.25) 
(4.26) 

(4.27) 
(4.28) 



Here we have used the definitions of M, N in (3.481,(3.49). Thus, the equations that we 



need to solve are simply (3.50 - 3.52) with the right hand side set to zero: 
1 



-9 



dr [V^9''''9''9xxM'] - qg'' {qM + wN) = 



1 



.dr [^/^g'^'N'] - wg^^ {qM + wN) = 
qg^^N' - wg^^M' = 0. 



(4.29) 

(4.30) 
(4.31) 



4.4 Gauge invariant variables 

As previously mentioned, it is desirable to work only with physical degrees of freedom, 
and this is accomplished by introducing gauge invariant variables which do not trans- 



form under (3.66). Within the context of the shear mode, one such gauge invariant 



combination is 

Z{r) = qM{r) + wN{r). (4.32) 

This can be seen by using the definition of the covariant derivative and the Christoffel 
symbols in Appendix [A} 

Vj^^Cv + V^o)^/. = + S.Cm - ^T^^M- (4.33) 



Using the presumed form of the vector ^ (3.67), we have 



(4.34) 
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and similarly 



(4.35) 



Thus, under the gauge transformation (3.66) the shear mode metric perturbations be- 
have as 



hta 



hta + iw^a 



or equivalently 



gxxN{r) 



9xxN{r) - iqia{r). 



(4.36) 
(4.37) 

(4.38) 
(4.39) 



It is clear that the variable Z{r) from (4.32) does not transform and is thus gauge 
invariant. 



It is quite straightforward to take combinations of the equations (4.29 - 4.31) and 



reduce them to a single equation for Z(r). Eliminating M from (4.31) results in 

wgxxZ'{r) 



N'{r) 



q^gtt + w'^Qxx 



(4.40) 



substituting this into ( |4.30 ) gives an equation which only involves Z. It can be written 

9rr{q^gtt + W^9xx) \ V^9xx ^,1 tt / 2 , 2 \ n 

Or — ; ^ - 9rrg [q gtt + W Qxx) Z = (4.41) 



-Z' 



XX J 

We now insert the incoming wave ansatz, 

Zir) = /(r)--/4^^(yo(r)+g2y2(r)+g^n(r)+0(g^)), (4.42) 
w{q) = W2q^ ^w^q^ ^O^q^) , (4.43) 



and expand (4.41) in powers of q. The Y functions must be regular at the horizon, and 



must also satisfy a Dirichlet boundary condition at r — > oo as explained in Sec. 3.2.5 



4.5 Solution for W2 and shear viscosity 

The lowest order equation for Yq is 



dr 



-99' y, 



0, 



(4.44) 



which is easily solved as 



lo(0 = ai 



fir') 



-dr' 
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(4.45) 



y^-gir')g^-{r') 

where ai is a constant, and the Dirichlet boundary condition has been applied at the 
boundary r ^ oo. Notice that the function Yq already meets the requirement of regu- 
larity at the horizon. Near the horizon, the integrand approaches a finite value due to 



(2.24), (2.25), and (3.57) and thus there are no divergences which need to be removed. 
For convenience, we define the quantity 



IW2 

4^' 



Then, expanding (4.41) to the next order in q gives an equation for Y2, 
dr 



-gg 

f 



Wo 



Yi + ^Y^-Cdr{Yolog[f]) 



-gg 



Performing a single integration gives 



^2 + y>^'-e5.(yoiog[/]) 



/ 



-gg' 



a2 



Yn 



-gg''dr' 



0. 



(4.46) 



(4.47) 



(4.48) 



Perhaps the simplest way to implement the requirement of regularity of Y2 is to 
require that, l^'(r — > tq) should be at most logarithmically divergent. To this end, let 
us examine (4.48) near the horizon. To leading order we have 

1 



Wn 



709^^ (?'o) 



Y^iro) - CYoiro) 



+ O(log(r-ro)) = 0. 



(4.49) 



In deriving the above expression, one can estimate the divergence of the right hand side 



of (4.48) by simply inserting the near horizon behavior of each function. For example, 

1 



the integrand 



Yo{r')y^^^g'\r') O 



(4.50) 



J.' — ro ^ 

so that after the integration over r', this term is at most logarithmically divergent, and 



does not contribute to the leading order in (4.49). By assumption, Y2 is regular at the 



horizon, thus the coefficient of the leading divergence in (4.49) must vanish. Recalling 



the definition of ^ from (4.46) we obtain the solution for W2'. 

-z7o5^''(ro) lo(''o) 



'W2 



AttT 



(4.51) 
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Using the definition of the temperature (2.391, and the solution for Yq (4.45), we find 

. \/-9{ro) 



W2 



ro 



fir') 
V-9{r')9"{r' 



-dr'. 



(4.52) 



Comparison with (|4.1|) and (|3.15|) gives a formula for the shear viscosity 



Dr 



JL 

Ts 



V-9{ro) 



\/-9{r')g''"''{r' 



-dr' , 



(4.53) 



which was already derived using similar methods in |115| . and from the membrane 
paradigm in [H 27] . Agreement with these results in a reassuring check on our method- 
ology. Explicit evaluation of this formula for any gravity dual which satisfies the equa- 



tions (4.29 - 4.31 ) results in 7]/s = l/47r, as should be expected given the general proofs 



of universality |l6l |l5] . 

It is worth reemphasizing that this formula has only been shown to be valid for 



gravity duals which satisfy (4.29) - (4.31), as well as the relevant background equations. 
There has been some confusion in the literature regarding this point. For example, 
one could try to use this formula to get the shear viscosity for a theory with a non- 



vanishing chemical potential (e.g. AdS Reissner-Nordstrom) . Application of (4.53) in 
this case results in r]/s ^ l/Air. The reason for this discrepancy is that such metrics 



are not generated by scalar fields alone, and hence the equations (4.29) - (4.31) are not 
satisfied. Of course, a more careful treatment of the perturbations in the case where a 
chemical potential is present gives the expected result, rj/s = I/Att \118\ 1119] I120| IT2T] . 



A second example where confusion may occur is application of the formula (4.53) 
to phenomenological models such as the soft-wall model. In the Einstein frame, these 



metrics are usually a warped version of anti-de Sitter space (cf. (4. 13), (4. 17)). Naively 



plugging such a metric into the formula (4.53) again results in r]/s ^ l/Air [T1 I122] . The 



reason for this strange result is exactly what was discussed in Sec. |4.2.2[ these metrics 
cannot be generated by scalar fields alone, and hence violate the assumptions inherent 



in the derivation of (4.53). 



In summary, if application of (4.53 ) to a background based on Einstein gravity results 
in a value different than l/47r, this result should not be viewed as a violation of the 
viscosity bound. Instead, it should be viewed as an inconsistency. All classical Einstein 
gravity duals saturate the conjectured viscosity bound as shown in 
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4.6 Solution for Tgi 



Liear 



We now proceed to the next order in q. By expanding the gauge invariant equation 



(4.41 ) to order (7^, one finds an equation for I4, 



Yi + wi^-f-idr[\og[f]Y2) 



59" (i2-eiog[/]ib) = o. 



+ ^0 ( ^ + ) - e ( ^ + ?1 5, (iog[/]yo) + fa, (iog[/]2yo; 



(4.54) 



Again, the requirement of regularity of I4 is best implemented by requiring that there 
are no power divergences in Yl near tq. Integrating once, we have 



n + -if -ea.(iog[/]y.) + yo'(^ + '"^"^^ 



r) 



^{j + '^)dr {log[f]Yo) + pr (log[/]^yo) 



/ 



-99' 



as 



-55" {Y2-Clog[f]Yo)dr' 



(4.55) 



We should now expand this equation near the horizon in powers of (r — ro). Before 
proceeding, it is useful to know the near horizon behavior of I2 and Yq. By construc- 



tion, lo(^o) and Y2{rQ) are finite. It is easy to verify from (4.45) that Yq^tq) is also 



finite. However, from (4.49), it is clear that 5^'(ro) could potentially be logarithmically 
divergent. Near the horizon, this equation takes the form 

Y^iro) - ^Y^iro) log(r - ro) + 0{1) = ^^^f^ ^ ^ 

The integrand on the right side behaves as 

\/-g{ro)YQ{rQ) 



pro 

/ ^g^'Yodr'. (4.56) 

J 00 



+ 0(1) 



(4.57) 



- ro) 

and is thus potentially logarithmically divergent upon integration. Taking this into 



account, (4.56) becomes 



F2'(ro) + log(r - ro) [-^^^'(ro) + 7r5""(r-o)lb(ro)] + 0{l) = 0. (4.58) 
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However, using (4.51), (4.46) and the expression for the temperature (2.39) we see that 



the term in brackets vanishes. Thus, we have shown that Y2{rQ) is also finite, along 
with Y2{ro),Yo{ro), and Y^{ro). 



With this knowledge, we return to (4.55) and examine its behavior near the horizon. 
For simplicity, we use the notation 



/(r ^ ro) « /o(r - ro) + /i(r - rof + 0{r - r^f 



(4.59) 



In all, (4.55) becomes 
yi^ro) + 7 ^2 



JO 



C^lo(ro) 

JO 



+ 



log(r - rp) 
r-ro 



fYoiro) - c'^Y^iro) 
Jo 



1 



r 



ro 
W2W4 



W2 
Jo 



Y^{ro) - ei2(ro) + ^ (Y^{ro) - 2^Y^{ro) 
Jo \ JO 



, Y^iro) - e^yo'(^o) - C—Yoiro) \ + 0(log(r - ro)) 

Jo Jo W2 



(4.60) 



The requirement of regularity of I4 at the horizon means that each of these divergent 



terms must vanish. Using (4.51), one finds that the leading and sub-leading divergent 
terms vanish. The coefficient of (r — vq)"^ does not vanish on its own; requiring that 



this term vanish we obtain an expression for W4. Again, employing (4.51), we find 

yiiro) Y^{ro)Y2{roy 



.W4 

i — 

W2 



^shear 



+ 



1 



Yoiro) YoiroY 



Yo{ro) 



47rTY^{ro) 



Yoiro) fo 



Y^iro) 



(4.61) 



This expression can be simplified somewhat by recalling a special property of our metric 



(4.12). Using the explicit expressions for Ft and given in the appendix, we see that 
our backgrounds have the property that yj—gg" oc ///'. Thus, 



yo'(r)«/'(r). 



Using this fact, and the expansion for / (4.59) we see that 



Y^'jro) 
Y^iro) 



2/1 

fo' 



(4.62) 



(4.63) 
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which ehminates the last two terms in the expression for Tshear given above. We can 
then write the expression for Tghear rather succinctly as 



Tshear = " ^ lim Or 



Y2ir) 

Yo{r) 



+ 



1 



AttT 



(4.64) 



All that remains is to get an explicit expression for Y2{r) in order to express Tghear in 
terms of the metric only. 



Returning to (4.48), one can solve for Y2/YQ: 



I2 _ 02 
Yq ai 



(4.65) 



It turns out to be useful to integrate the term involving the nested integrals by parts. 
With this simplification, we have 



YZ "2 , ^1 r„i 1 

Yo ai ai 



-gg^'Yodr' 



+ 



1 



00 J 



Taking the derivative. 



dr 



Y2 

Yo 



fYo + Yi 



I r-r /vi„,,2 



Yo^ 
f 



ai 



dr' , 



(4.66) 



(4.67) 



and using the definitions of Yq (4.45), and W2 (4.51) this becomes 
dr 



Y2 
Yo 



fYo 



+ ai 



Y^jr) fifoYoirol 
yo(r)2 V4vrryo'(ro) 



9r 



-9 



1 + 



Yo^ 



V-9tt9rraiW2 



(4.68) 
dr'. 



Defining the quantity 
D{r) 



1 V-9{r) 



V-9{r) 



fir') 



^ ^ =Yo{r) = ^ ^ = / ' dr', (4.69) 

«i V-9tt{r)grr{r) \/ -9tt{r)grr{r) Jr \J-g{r')g"{r') 

and noting that 



W2 = -iD{ro), 



(4.70) 
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we have 



dr 



'Y2 



fYo 



9r 



1 



D(r 



i\2 



dr. (4.71) 



-g V D{rQ)\ 

The first two terms of this expression can be simplified by again invoking the special 
property of our background that Ft{r) = Fx{r) which means that ^J—gg^^ oc /'// so that 
(4.45) gives YQ{r) oc 1 — f{r). Using this fact, and the universal fact that W2 = —i/^TiT 

1, we have 



for all Einstein gravity duals 



1 



fir) 
4ttT J /(r)-l' 



(4.72) 



Then, with the use of (4.59), we have 

/o 



lim dr 

r^ro 



Y2 

Yo 



+ ai 



(4^r)2 
1 



A. 

yo'(ro) V4vrT 



i\2 



Inserting this into (4.64) gives 



1 



''"shear '^1 



/o 



yo'(ro) 47rr 



^•0 



i\2 



W) 

D{ro) 



D{roy 



dr' . 



dr. (4.73) 



(4.74) 



As a final step, insert the definition of Yq (4.45) and substitute in for the temperature 
using (2.39). This precisely determines the quantity Tshear which thus completes the 
shear dispersion relation up to order q^, 



''"shear 



1 



D{r) 
Wo) 



(4.75) 



where D{r) is defined as in (4.69). 



The formula given in Eq. (4.75) is the main result of this chapter. This is the 
same formula that was derived in [1] within the context of the black hole membrane 
paradigm. Here, we have shown that the formula can also be derived using AdS/CFT 
and the method of gauge invariant equations. This further illustrates the agreement 
between these two methods, at least in the case of zero chemical potential. For a more 
in depth discussion of this agreement (which is not fully understood at present) see 
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4.7 Applications 

Let us consider a few simple applications of this formula. All metrics generated by r 
dependent scalar fields will have the aforementioned property that Ft = Fx, and thus 
^—gg^^ (X f/f. In addition to this "automatic" constraint, let us examine metrics 
which satisfy 

Here, C2 is a constant. We also assume that the metric is asymptotically AdS so that 
/(r) 1 as r — > oo. We will explain the motivation behind this constraint in subsequent 
chapters; for the moment, imposing the constraint simply allows us to get analytical 



solutions for a special class of metrics. Using (4.75) for this special class of metrics we 
have 

rshear = ^Hr. (2 - ^) , (4.77) 

where Hn{a) is the "Harmonic Number" defined as 

Hn{a)= / dx. (4.78) 

Jo ^ - X 

The results for r//s and Tghear for this special class of metrics is represented graphically 



in Fig. 4.1 



Two notable special metrics are the Schwarzs child anti-de Sitter black hole, and the 
Dp-brane metric. The first of these is a black brane type metric which is a solution to 
Einstein's equations with a negative cosmological constant. This (p + 2 dimensional) 
metric can be written 



(v>2 T 2 

SA,S = j2[-mdt' + dx,dx^]+^^ (4.79) 



/(r) = (4.80) 



r 



Here, L is the usual AdS curvature radius. Clearly, this fits the constraint (4.76) with 
C2 = {p + l)/2, thus for this metric we have 

^i^-i^^.(^). (4.81) 

As a second example, the Dp-brane metric arises from a stack of N Dp-branes in string 
theory with ^ 1. The metric exists in 10 dimensions, but can be reduced to an 
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effective (p+2 dimensional) metric using dimensional reduction |47| 1101] . The effective 
metric is a solution to Einstein's equations where the metric is supported by a single 
scalar field, and the scalar potential is a single exponential 



U{<t>) 



{7-p){9-p) 
2L2 



(9-p) 



The metric is: 



9-p 
V 



[-f{r)dt^ + dxjdx^] + 



L 



-8p+9 



dr^ 

W) 



r 



(4.82) 

(4.83) 
(4.84) 



This fits our special class of metrics with C2 = {7p — p^) /(9 — p). (Note that p < 7, or 



else the metric is not asymptotically anti-de Sitter.) Applying the formula (4.75) we 
have 



^DP 
' shear 



1 



4ttT 



7- 



P 



(4.85) 



The results (4.81 ) and (4.85 ) were first presented by Natsuume in jl05j . A few comments 



are in order about these applications. First, before the formula (4.75) was presented 



in [T], the quantity Tghear had been computed in several special cases, most notably in 



|89] |90] . The applications above (4.81) and (4.85) reproduce all known results in the 
literature, which is another check on our calculations. In addition, the formulas above 
give new information. A particularly curious feature is that t^^^ = for p = 5, and 
"^shear ^ p > 5. It has been known that the Dp-brane metric exhibits unusual 
features for p > 4 in the sound channel. (For example, for p = 5, the speed of sound 
vanishes, while for p > 5 the speed of sound becomes imaginary |101j ). The analysis 
given here shows that these features are also present for the shear mode. 

Furthermore, once we have completed the analysis of the sound mode, the compar- 



ison with the formula (4.77) gives us some information about the third order hydrody- 



namic contribution to the shear mode as mentioned in Sec. 3.1.3 This will be explained 
further in the next chapter. 

Finally, perhaps the most obvious feature of these applications is that the quantity 
T"shear is not Universal in the same way as the shear viscosity; application of the formula 
to different metrics yields different results. 
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4.8 Summary 



In this chapter, we have apphed the prescription which was outUned in Sec. 3.2 to 
shear perturbations of a black brane background supported by scalar fields. We have 
thus been able to determine the shear mode dispersion relation 'w{q) up to order q'^. 



This dispersion relation is parametrized in (4.1) with two unknown quantities which 
we denote Dfj and Tghear- The formulas (4.69) and (4.75) are the main results of this 



chapter, and they allow one to compute these unknown quantities if the background 
metric is given. 

The formula for Z)^, which was first derived in jJT], allows one to compute the shear 
viscosity, since it is the only hydrodynamic transport coefficient which appears to lowest 
order in the shear dispersion relation. Application of the formula for D^i to any Einstein 
gravity dual (in the limit of infinite coupling) yields rj/s = l/47r, and is one of the main 
inspirations for the conjectured shear viscosity bound. The quantity Tghear does not 
display the same universal behavior. 

It is important to make a distinction between the dispersion relation and the hydro- 
dynamic transport coefficients which would appear in the derivative expansion of the 
energy momentum tensor. To lowest order, the shear dispersion relation contains only 
the shear viscosity, but at the next order (g^), the quantity which we call Tghear contains 
more than one transport coefficient. (In fact, it contains contributions from third order 



hydrodynamics as explained in section 3.1.3). At this point, all we can state is that the 



combination of transport coefficients which comprise Tghear is not universal in the same 
way as the shear viscosity. If and when a consistent theory of third order hydrodynam- 
ics is formulated, one could revisit these results to see whether additional information 



about these new transport coefficients can be gained from the formula (4.75). 



57 



r]/s 

Tshear(4 71 T) 



p=2 p=3 




Figure 4.1: Plots of the dimensionless quantities: rj/s and 47rTrshean ^-s a function 
of the free parameter C2 appearing in the metric (4.76). Each plot corresponds to a 
different number of spatial dimensions p. This range of values of C2 was chosen because 
outside this particular range, unphysical behavior is found in the sound mode. This will 
be explained in later chapters (cf. (5.861,(5.87)). The metric is conformal at the upper 
end of the range of C2 when C2 = (p + l)/2. Note the universality of the shear viscosity 
to entropy density ratio is evident, while the value of Tghear varies with both C2 and p. 



Chapter 5 



Sound Mode 



5.1 Introduction 

In this chapter, we will perform an analysis similar to that of the previous chapter 
within the context of a different hydrodynamic mode, the sound mode. In general, the 
prescription detailed in Sec. |3.2| will be applied in exactly the same way, but analysis is 
more complicated here due to the fact that the matter perturbations are coupled with 
the metric perturbations. 

In [37], a formula for the shear viscosity to entropy density ratio ry/s was derived 
which is applicable to a wide variety of gravity duals. Application of this formula to 
known Einstein gravity duals always resulted in rj/s = l/Air (in the large N/infinite 
t'Hooft coupling limit), and eventually led to the celebrated viscosity bound conjecture 
that r]/s > 1/4tt for all physical substances. Clearly, a similar general formula for the 
bulk viscosity would also be desirable in the hopes that other universal behavior might 
be discovered. Such a formula would also be useful when comparing the predictions of 
a phenomenological gravity dual to the lattice data |123| I124j . 

The work detailed in this chapter is partially motivated by the goal of obtaining 
such a formula, though such a formula has not yet been developed. This work should 
be viewed as a first step towards this goal. Here, we consider sound mode fluctuations 
of a dual gravity theory supported by a single scalar field, and then generalize to the 
case of multiple scalar fields. We attempt to be as general as possible by not specifying 
the background profiles of the scalar field(s) except in certain special cases. 
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Recent works |37| 1125] also explore the sound mode in phenomenologically motivated 
single scalar models. The results presented here are complementary to these papers, 
though we allow for the possibility of multiple scalar fields, and also do not restrict 
ourselves to five dimensions. 

The main results of this chapter are twofold. First, we derive a set of sound mode 
gauge invariant equations which are applicable to any Einstein gravity dual which is 
supported by r dependent scalar fields. These equations are, in general, complicated 
but could be solved numerically to determine the sound mode dispersion relation for 
certain phenomenological models of the quark-gluon plasma. Second, we solve these 
equations analytically for the sound mode dispersion relation (up to order q^) for a 
certain special class of metrics. 



5.2 Background field and potential 



We again work with the same set up as in the previous chapter. We will examine sound 



mode perturbations of a general black brane metric (2.23) with presumed near horizon 



behavior (2.24 - 2.26). We assume that the matter supporting the metric is scalar in 



nature, so that the action, background equations, and energy momentum tensor are 



given respectively by (4.2), (4.4), and (4.5) 



Because the sound mode perturbations transform like scalars under rotations (see 



Sec. 3.1.2), the perturbations of the scalar fields will be coupled to the metric pertur- 



bations. For this reason, the sound mode requires more information about the matter 
perturbations than the shear mode. In particular, we will need information about the 
scalar potential. 

The background equations of motion give a relationship between the scalar field 
profile and the potential. 



dUiUr)) 



--dr [^/=5/>:(r)] . 



(5.1) 



As usual, the prime denotes derivative with respect to r, and D^^'^ = V^'^^Vj^Qj. The 
covariant derivatives are defined in Appendix [A} 



5.3 Perturbations and linearized equations 
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We now introduce fluctuations of tlie fields on this background g^u gfj,u + /i/^i/ and 
(j) ^ (f) + 5(f), and assume the usual time dependence 



64>{t,z,r) 



J{qz-wt 



(5.2) 
(5.3) 



The linearized Einstein equations for the sound mode and an arbitrary set of matter 



fields were derived in Sec. 3.2.3 In order to apply these equations to the case at hand, 
we must evaluate the perturbed energy momentum tensor for the case of the scalar 
fields. The energy momentum tensor was expanded to first order in the perturbation 



in (4.20). For the sound mode, we require the following components of the energy 
momentum tensor: jT*^"" ,T^^ ,T*^ ,T^^ ,T^'^ ,T'^'^ . As in previous chapters, we use the 
subscripts (0) and (1) to distinguish background quantities from those that are 0{hfj_u) 
or 0{6(j)). A summation over the repeated index k, which labels the particular scalar 
field, is implied. 



First, let us examine with m,n ^ r. Using (4.20), we have 



rpmn 
OVrLrp-|.2-t(l) 



Noting that 



16ttGp+2T, 



(0) 



„mn f 

-g J-d,, 



and recalling the definitions of A,B,C, and D ( 3.53 )-( 3.56) we have 

167rGp+2(r(*^) + r**)A) = -g''6C^, 

= -g^^6£., 



167rGp+2 \ T^{) + T^Q)C 
16^Gp+2 (r*f) + Tfop 
The right hand side is straightforward to evaluate 

5Crh - 



0. 



(5.4) 
(5.5) 

(5.6) 

(5.7) 

(5.8) 
(5.9) 



dU 

k 



+ 5U 



9 



k)'4)'h + 



(5.10) 
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Secondly, examine TJ^^: 
This gives the relations 

Finally, the only remaining component to evaluate is T^{y evaluating (4.20) gives 

[2 {g'''fdriS(l)k)dr<Pk-g'''SC^ 



which simplifies to 



167rGp+2T[[) = {gnHHky^'k-9 



\' J,' „rr 



, dU 



(5.11) 

(5.12) 
(5.13) 

(5.14) 
(5.15) 



The insertion of these relations for the energy momentum tensor into the general sound 



mode equations (3.59 - 3.65) results in the following set of seven linearized Einstein 



equations. Here we have inserted a line above each equation that allows the reader to 
see which Einstein equation is the source of the following differential equation. 



til) 



fdr 



-99 



G 



-{(p-l)B' + C') 



dU 

+ 9rrTr-{0(Pk) 



(5.16) 



+ </.',(<5c/.fc)'-(p-l)(zVr5""i? = 0, 

x;Gf)=g(-8.G,,2T:« 



i=l 



1=1 



dr 



^g^^A'] + -^dr[V^gg^'{{p-2)B' + C')]+gr 



dU 



-9 

7 

9rr9 


G 



w 



{A + {p- 2)B) - —{{p- 2)B + G) + 2qwD 



f 



(5.17) 



zil) 



dr 



dU 



-I V-9 ocpk 

- {p-l)w^grr9''B = 0, 



(5.18) 
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^(1) 



-g{a 



xx\p+l 



-dr 



G 



drifD) 



qrr{Ft - F^)D + qw{p - l)5rr5"^ = 0, 



G 



r(l) 



dU 



(5.19) 



(5.20) 



2grr{w'g'\ip-l)B + C) + q'g''''{A+ip-l)B)+2qwg''''D) = 0, 



Gr^ — — 8TrGp-\-2Tr^ 



w 



yffdr 



Vf 



{{p-l)B + C) 



qfD' + w(t)'i. 



(5.21) 



G 



zil) 



-8TrGp+2T' 



z(l) 



Vf 



dr 



w 



fA + jdr [fD] + q{p - 1)B' + #U'^</'fc) = 0. 



(5.22) 



In addition to these equations, one must also examine perturbations of the background 



scalar equations (5.1). Using the formulas given in the Appendix, especially (A. 19) the 
equation for each scalar field (in the radial gauge, as usual) is 



9r 



-dr [^/=5/'^(50.)'] + h^[H' - g„ {w^g'' + q^g'''^) {5<l^i) = 9r 



-g ^ - - ^ - _ _ , , , _ d(j)id(l)k 

Again, a summation over k is implied and we have defined 



i5<Pk)- {5.23) 



H{r) = A{r) + (p - l)B{r) + C(r). 



(5.24) 
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5.4 Gauge Invariant Equations 

There is still residual gauge freedom in these equations under the infinitesimal diffeomor- 



phism (3.66). The scalar field perturbations also transform under a gauge transforma- 
tion. Such a gauge transformation is essentially the freedom to choose ones coordinates. 
Imagine a small coordinate shift 

xf" ^xf"- (5.25) 

Then, the background scalar field will be altered 

(^(x^) ^ (t>{x'') - ^f'df,(l){x''). (5.26) 

The latter term appears to be a perturbation, but it is not a physical perturbation; it 
only appears due to a change in coordinates. Thus, under such a transformation, the 
scalar field perturbations behave as 

64>i ^ - i^'d^ct^i = Hi - ^'(t>-- (5.27) 

Let us now examine how the different sound mode perturbations transform under this 
gauge transformation. Recall 

V(f)6 + V(°)Cm = d^Cu + d,^^ - 2r^,^A, (5.28) 
the Christoffel symbols given in Appendix [A] and the presumed form of the vector ^ 



(3.67). Then, it is straightforward to evaluate 

htt ^ htt-2dtit + 2Tiix = htt + 2iwit-g'''g'uir, (5.29) 

A ^ A + liwo^'it-g'^^VLlgu^. (5.30) 

In a similar manner, 

B ^ (5.31) 

C ^ C-2ig5""e.-/'^^L[fc]^r-, (5.32) 

D ^ D + iwg''i,-iqg''it. (5.33) 



Given the transformation rules (5.27), (5.30), and (5.33), it is easy to verify that 
the following gauge invariant combinations first presented in |101] transform only into 
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themselves under such a diffeomorphism (i.e. Zq — > Zq and Z^ Z^). The choice of 
these gauge invariant variables is not unique. 

2 9u{r) 



Zo{r) 



-f{r) [q^A{r) + 2qwD{r)\ + w'^C{r) - 

<t>'i{r) 



Z^i{r) = 5(j)i{r) 



T^L [gxx{r)] 



B{r) 



+ w 



B{r) (5.34) 



(5.35) 



By taking an appropriate combination of the equations (5.16 5.23), one can arrive at 



coupled second order equations for the gauge invariant variables. Because of the com- 
plexity of the sound mode equations and the gauge invariant variables, reducing the 
linearized perturbation equations down to a set of gauge invariant equations is a very 
difficult task. The equations presented below were arrived at by "brute force" , by elim- 
inating the variables A, C, D, 5(j) one at a time, and then by showing that the remaining 
coefficient of B vanishes by the background equations of motion. 



5.4.1 Equation for Zq 

First, define the quantity 



a(r) = \{p-l) + 



T^L[gtt{r)] \ pw 



VL[gxx{r)]) f{r)- 
Consider the following combination of the Einstein equations ( |5.16 - 5.23): 



(5.36) 



/ 1^ [(|5l6| - ^J7\ + (|5l8( + ^(|K20)] + [jsle) + ^(|K20)] (^y^ - 2qw{5.19] 



(/'[^(pof + ilslsl)] \ + 2 {qfVL [a] ([532|) - wVl [af] } . 



(5.37) 

For simplicity, for the moment we assume there is only one scalar field. After a long 
calculation, one can show that this reduces to 



Z 



1 -^yrr 



{Ft - F,) {w\B -C) + Zo) - Bq^fAi = 0, 



(5.38) 



where 



9rr 2 f2 



-9a 



-z' 



2Z^4)'f ad, 



"1 /w^ 



+ Zo {VL[f]VL[fa] - g„ {w^9'' + q'^"")) + 
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and 



9r 



1 



{Ft - F,) 



(2grr{Ft — Fr 



l\2 



(5.40) 



Here, we have been able to eliminate the potential using (5.1 ), since it only appears in 
these equations via dU/d(j). Clearly, Ai vanishes by the background equations of motion 
(4.12) and ( 4.10| ), leaving a differential equation involving only the gauge invariant 
variables, namely Zi = 0, 



vr 

99 y, 



ay 



+ Zo(PL[/]PL[/a]-<7rr(W* + gV")) + (5.41) 



2Z^07 ad, 



1 I w 



a 



f 



+ 



pT>l[9xx] 



-99 



TT i/1 



In this special case of one scalar field, one can write these equations only in terms of 



the metric components because of the relation (4.10). 

The generalization to multiple scalar fields is straightforward. Since all of the equa- 



tions used to derive (5.41 ) depend only on the sum of all scalar fields, it is clear that all 
we need to do is include a summation on the last term. An explicit calculation (using 
the same combination of background Einstein equations) reveals this is correct. The 
generalization is 



'-9 



VT 

99 ry, 



k=\ 



+ Zo {piX!\ViXfa\ - g„ [w'g'' + q^g^^)) + 
1 (n? AI q^ViXn 



+ 



V^l\9x 



~99"<^k\ 



0(5.42) 



5.4.2 Equations for Za 



We now proceed to derive the other gauge invariant equations. There will be an addi- 
tional gauge invariant equation for each background scalar field present. The relevant 
combination of equations is: 



( |5.23[ ) + 

+ 



vDl\9'. 
2 



(5.17) - (5.18) - (5.16) - -(5.20) 



a 



-a. 



T^l{9x 



(5.21) -g(5.22) 



(5.43) 



Another lengthy calculation reduces this combination to 
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Z2 - BA2 = 0, 



(5.44) 



where 



Z2 ^ -^dr[V^rZ'^,]-Z^,grr{w'g''+q'g--)-grrY.Z, 



-9 

20'. 



(f>k] 



k=l 



p'DL[gxx\a 
2 



oi\ff 



dr 



E 

.fc=i 



Zo 



'DlIqx 



w 



5, 



and 



Ao = 



+ 



xx\ 
I 



E 



92 [/ . 



(5.45) 



k-9r 



PVlIQxxY 



9rr 



dr 



+ 



+ 



a'DL[gxx_ 



■dr 



F>L[gx 



2grr ( 2q 



T^l[9xx] V « 



-dr 



+ c^'A{Ft-F,). 



(5.46) 



Again, A2 vanishes by the background equations of motion (4.10), (4.12), and (5.1). 
Thus, we are left with an equation only containing the gauge invariant variables, Z2 = 0: 



gr 



-g 

2cp' 



dr [V^g"'Z'^,] - Z^,grr {w^g'' + g'''^) - grr Z 



<j>k- 



k=l 



9 

' -7 



+ 



dr 



T^dgx 



dr 



z^ 



0. 



(5.47) 



When only one scalar field is present, (5.47) simplifies. This equation contains the 
combination d'^U/d(j)id(j)k- In the case of one scalar field, one can use 



(fU{(j){r)) 1 



(/)'(r) 



dr 



dU{(j){r)) 



(5.48) 
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which, along with the expression for U'{(j)) (5.1), allows us to eliminate the potential 
completely in favor of the field (j). 



\ 1 













-g '- a ^ „ , ^ 



w 

T 



{i^Ligxx]) 



-gg 



(5.49) 




The equations (5.42) and (5.47) are the most general gauge invariant equations for 



any black brane background supported by r dependent scalar fields. If only one scalar 



field is present, one can use (5.41) and (5.49). To determine the dispersion relation 



one needs to solve the above equations perturbatively in w,q as explained in previous 
chapters. Some terms can be neglected when considering these equations to 0{w,q), 
but one should keep in mind that Zq is 0{q'^), Za^ is 0(1), and a is 0{q^). 



5.5 Special case solution 
5.5.1 Solution for 

Even for the case of a single scalar field, the gauge invariant equations are quite com- 
plicated and we have been unable to find a general analytic solution. However, analytic 
solutions do exist for certain special backgrounds, as we now show. 

Consider the special case of one scalar field with profile (j){r) = i^\og[gxx{f)] + </'o 



where k and 0o are constants. In this special case, the equation (5.49) reduces to 

grr 



-g 



dr [V^gg^Z'^] - Z^grr {q^ g^"" + w^g'') = 0. 



This can be more easily seen by re-writing 

/0' 



a, 



-gg 



Vl W^gg"VL\gxx\\ 
2grr{Ft 



(^L[fc]) 
/ 

T^dgxx] 

0, 



Fx 



fdr 



p{VL[gxx]f 

(/)'(r)2 " 

viPL[gxx]f . 



(5.50) 

(5.51) 

(5.52) 

(5.53) 

(5.54) 
(5.55) 
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where we have used (4.10) and (4.12) in the last two steps 



We now proceed according to the prescription outlined in Chapter |3] by applying the 
incoming wave ansatz 



^o(r) 
w{q) 



f(r)-^ (Yoir)+qYi{r) + ...) 
f{r)-^{Y^o{r)+qY^i{r) + .. 
wiq + W2q'^ + W3q^ + ... 



(5.56) 
(5.57) 
(5.58) 



with the condition that all Y functions are regular at the horizon. Inserting this ansatz 



into (5.50), and expanding the result in powers of q, we have the following equation 
which is good up to O(q^): 



drlV^g'^^.j+qdr 



-99 



Y' 



2^ 



(5.59) 



It should be noted that we have omitted some terms which are proportional to {Rq — R^) 
since they vanish by the background equations of motion. Solving this order by order 
in q is now quite simple. The solution for Y^q can be written in terms of an integral 

9rr{r') 



Y^oir) = Co + ci 



\/=5F) 



dr', 



(5.60) 



but this integral is logarithmically divergent at the horizon by (2.25). Thus, the as- 
sumption of regularity on the Y functions leads to Y^q = cq. Finally, this constant must 
be set to zero by the Dirichlet boundary condition at infinity. Plugging l^o = into 
the next order equation, one also finds l^i = 0. Proceeding now to higher orders in q, 
one finds that the equations always reduce to the same as that for Yq, and as a result 
= to all orders in q. With this simplification, the gauge invariant equation for Zq 
that must be solved is 



9r 



-99 



2 f2 ^0 



+ Zo{VL[f]VL[fa]-grr{w^9'' + q'9'''')} = 



(5.61) 



5.5.2 Constraints on the metric and scalar potential 

Next, one must solve the equation for Zq using these boundary conditions with the 
knowledge that Z(f, = 0. Before doing so, it is useful to pause and ask what type of 
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metrics these results will be applicable to. There are three unknown metric functions, 



gtt,gxx, and grr, but we have two constraints on them, namely (4.10) and (4.12). These 



two equations allow one to solve for both gu and g^ in terms of gxx- Solving (4.12) 
gives grr in terms of the other metric components and one integration constant ci, 

grr{r) = cigxx{rr+'f{r)VL[f{r)f. (5.62) 

For convenience, we define the constant of proportionality between (p{r) and log[g(xx] as 



(r) 



where (po is an integration constant. 



Klog[gxx{r)] + 4>o 
±,r-z{p+l-2c2). 



(5.63) 
(5.64) 



fir) 



C3 



1 — C4 



(5.65) 



Taking this solution, plugging into (4.10) and solving for gu gives /(r) in terms of 
gxx and three constants C2, C3 and C4. These constants could in principle depend on ro, 
but are independent of r, 

1 y^' 

^gxx{r)J _ 

Imposing the requirement that the function / must vanish at the horizon determines C4. 
Furthermore let us assume that the metric is asymptotically AdS, so that for r — > 00, 
gxx{r) oc r". Where n is some positive exponent. Then we must have C3 = 1, and C2 > 
for the correct asymptotic behavior of /. This determines / up to a single constant, 

' 9xx{ro)^ 



1 



9xx{r) 



(5.66) 



which is exactly the same metric constraint we considered previously (4.76). In the 
previous chapter, we simply imposed this constraint by hand. Here we have shown that 
the constraint arises naturally by first imposing a particular condition on the scalar 
field which simplifies the gauge invariant equations. The Hawking temperature for this 



metric can be found from (2.39) as 



1 



cigxx{roY 



(5.67) 



(47rr)2 

It is also interesting to ask what sort of scalar potential can generate such a metric. 
One can determine the form of the potential by considering the following combination 



of background Einstein equations 

An explicit computation of the left hand side for our special metric gives 



[/(</.(r)) 



P 



2ci C2 
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(5.68) 



(5.69) 



-2C2 



exp 



(5.70) 



Using (5.63), we find 

P 

= gxx{ro) _ , 

2ciC2 V p 

Thus, the sort of metrics we are considering are those generated by a potential which 

contains a single exponential, similar to the Chamblin-Reall backgrounds |126| examined 

in |36| I125j . It is amusing that we have arrived at the same metric considered in these 

references, as our only guiding principle has been to search for a class of metrics which 



simplify the gauge invariant equations (5.41 ), (5.49 ). 



Finally, note that since the potential must be independent of temperature, it is 
required that the constant C2, and the combination cigxx{'''o)'^'^^ must themselves be 
independent of tq. 

5.5.3 First order hydrodynamics 



One must now go back to (5.41), insert Z^j, = and the incoming wave ansatz (5.58) 



and expand the resulting equation in powers of q. In doing so, one should take into 



account the constraints (5.62) and (5.66). The first of these constraints allows one to 
eliminate \/—gg^^ in favor of Dilf]. Once this is done, the lowest order equation for Yq 
can be written 



{Co-fip-c2)y 

{P-C2)f 



dr 



1 



—V' 



Co-f{p-C2) 



Yo = 0, (5.71) 



VL[f]iCo-f{p-C2)Y 

where 

Co = pwl-C2. (5.72) 
The general solution to that equation contains two arbitrary constants ki and ^2, and 



can be written 

Yo{r) = [pwl-C2 + f{r){p-C2)] 



ko + ki 



pwf 



C2-{p-C2)f{T'fV f{r') 



pwf -C2 + {p- C2)/(r')^y f{r') 



dr' 



(5.73) 
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as can be found by first making the ansatz Yq^t) = bo + 6i/(r), and then using the 
technique of reduction of order once this solution is found (cf. Appendix p|. The 



integral in (5.73) is logarithmically divergent near the horizon, and thus the requirement 
of regularity leads to ki = 0. Finally, applying the Dirichlet boundary condition at 
r — > oo, leads to 

1, (5.74) 



2C2 

p 

where we have assumed that /(r oo) = 1. 

Proceeding now to the next order in q, and substituting the solutions for wi and Yq, 
one finds the following differential equation for Yi: 



dr 



VL[f]{i + fy 



+ 



f 



(1 + /)= 



0. 



(5.75) 



The solution to the homogeneous part can be found using the same techniques as for 
the solution for Yq listed above. A particular solution to the inhomogeneous equation 
is obviously Yi(r) = Constant. This leads to the general solution, 

(/(/) + 1)2 /'(/) 



Yiir) = {f{r) -I){k2 + k3 



{f{r') - 1)2 f{r') 



dr' I + wiko 



i{p - C2) 



vrT 



+ 2pW2 

(5.76) 

The integral above is again logarithmically divergent, leading to the requirement that 
^3 = 0. Applying now the Dirichlet boundary condition at r ^ cxd leads to 



wiko 



i{p - C2) 



+ 2pw2 = 0, 



W2 



i{c2 - p) 



(5.77) 



(5.78) 



2ttTp 

Note that we require C2 < p to avoid exponential growth of the perturbations. Let us 
summarize the central results of this section: 



Yo{r) 
Yi{r) 

Wl 
W2 



yo {f{r) - 1) , 

yi (/(r) - 1) , 



2C2 _ 
p 

i{p - C2) 
2'kTp 



1, 



(5.79) 
(5.80) 

(5.81) 
(5.82) 
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Here i/q = A;o(c2 — p) is a constant, though its value is unimportant as it does not enter 
into any formulas for physical quantities, and we have relabeled k2 as yi for aesthetic 
reasons. 

Comparing these results to the dispersion relation expected from first order hydro- 



dynamics (3.27) 



W2 



e + P 



P 



p 2ri 



(5.83) 
(5.84) 



one gains knowledge of Vg (speed of sound) , r] (shear viscosity) , and C, (bulk viscosity) : 



ry/s 



1/4^, 

p 
1 



(5.85) 
(5.86) 

(5.87) 



Here we remind the reader that s = (e + P) /T is the entropy density, e is the equilibrium 
energy density, and P is the equilibrium pressure. Note that the conjectured bulk 
viscosity bound of Buchel |127] is saturated for metrics of this type. 



5.5.4 Digression: comparison with black hole thermodynamics 

It is interesting to note that the speed of sound can be derived solely from the metric by 
using relations from black hole thermodynamics. Such a formula was derived in Chapter 
[2] (cf. 2.45). Applying it to the special metric and using (2.39) and (5.67) results in 



p + 



d 
dro 



Log[ci] 



(5.88) 



^^05([fc(ro)]_ 

Recall that ci is independent of r, but could in principle depend on tq. At first sight 
this seems quite puzzling due to the presence of the constant ci which is not present in 



the formula (5.86). These two calculations of the speed of sound must agree. 



The resolution of this apparent contradiction lies in the observation made at the end 



of Sec. (5.5.2). In order to have a sensible potential which is temperature independent, 
both C2 and ciQxxi'ro)'^^^ must be independent of temperature (and thus also independent 
of ro). As a result, we have 



ci oc gxx{ro) 



-2C2 



(5.89) 
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where the constant of proportionahty is independent of r and rg. Substituting this into 



(5.88) provides 



2C2 

p 



1, 



(5.90) 



which is exactly the same as that computed from the gravitational perturbations. It 
is quite reassuring that these two calculations agree. It is interesting that this latter 
method is much more efficient for computing the speed of sound; unfortunately, it does 
not give access to the bulk viscosity and any higher order transport coefficients. 

5.5.5 Second order hydrodynamics 

It is the purpose of this section to extend the above calculation to the next hydrodynamic 
order, and thus to determine the next coefficient in the dispersion relation w^. 



Equation for Y2 

We have already shown that vanishes to all orders in q, so then it remains to return to 



(5.61 ), insert the incoming wave ansatz, expand in powers of q, and insert the solutions 



(5.79 - 5.82 ). Completing these steps, one finds the following differential equation which 



must be solved for Y2: 



V 



f 



Y2 + 



+ 

yo 



(5.91) 



xo + 



1-f 



gxx[ro)P 



where 



xo = 2[5wt-l 



1 



(5.92) 



In writing the above expression, we have replaced all occurrences of the constant C2 



which appears in the metric with wi due to the relation (5.81), and have removed the 



constant ci in favor of T using (5.67). 



Solution for Y2 



The associated homogeneous equation for (5.92) is the same as the homogeneous part of 



the equation for Yi, (5.75 ). The solution can easily be found by ffist making the ansatz 
Y2{r) = y2a(l — f{f))i s-iid then using the technique of reduction of order once this 



74 



solution is found. The general solution to the homogeneous part of the above equation 
is 

Y2h{r) = y2a(l - f{r)) + ^26 [(1 - f{r)) log[/(r)] + 4] . (5.93) 

Here the subscript h stands for homogeneous, and y2a, y2b are arbitrary constants. 

In order to find the general solution for Y2, we must now find a particular solution 
to the inhomogeneous equation. Since the homogeneous solution is known, one can 
construct a particular solution using the method of variation of parameters. For com- 



pleteness, this method is outlined in Appendix B.2 Using the notation from Appendix 



B.2 (and changing independent variables from x to r), we have 

yi(r) = I- fir), 
h{r) = log[/(r)]+ ^ 



1 - fir) 

ifr ^0 
/(I + /) (47rr)2 



, 2 , 9x.ir)^ I, 



/ 



gxxiroY 



After some simplification, one can write the particular solution as. 



^2p(r) 
Q{r) 



yo 



iAirT)' 



r(l-/(0) 



l + fir)Y fir) 



(l + /(z))3^^"^ 



xo + 



I -fir) J fir) 



Qir) dr 



2 , 9xxiz)^ ( f. 



9xxiro, 



p 



(5.94) 
(5.95) 

(5.96) 



(5.97) 

dz. 
(5.98) 



The term involving xq can be integrated directly, by changing variables from z to f. 
This simplifies the solution to 



Y2pir) 



Qir) 



xoyo 



+ 



yo 



r(/(0 



(4^r)2 (47rr)2 



l-/(r)\ m 



l + fif)J /(f) 



1) 



wf + 



l + fir)Y fir) 



I -fir) J fir) 

9xxif)^ 



9xxiroy 



Qir) dr 



dr. 



(5.99) 
(5.100) 



and so, the general solution for I2 is given by 

Y2ir)=Y2hir)+Y2pir). (5.101) 
To proceed, it is convenient to change coordinates. We define the coordinate u by 

'9xxiro)Y' /fc(ro)^'^"'+'^ 



u 



9xxir) 



9x 



(5.102) 
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so that 

f{u) = l-u'^, (5.103) 

and the horizon is now located at u = 1, and the boundary is located at u = oj^ In 
terms of the new coordinates, I2 is given by 

y2h{u) = y2au'^ + y2b u"^ log[l - U^] + 4 



Y2p{u) 



yo 



(2 - u') 



2\2 ru V 



w 



2 ^ ^~4/(l+«,2) 



V — Wi 



(5.104) 



dv du 



(47rr)2l J V?{l-u'^)J (2-^;2)2(l_^2) 

(5.105)' 

The inner integral in this latter equation can be written in terms of hypergeometric 
functions, but the result is not particularly enlightening, so we do not reproduce it 
here. The full analytical form of this function is not needed to determine the dispersion 
relation. 



Boundary conditions 

The first boundary condition which must be applied on Y2 is the condition of regularity 
at the horizon. In order to do so, one must extract the coefficient of the logarithmic 
divergence in the particular solution Y2p. To do so, we first expand the integrand in 
powers of {u — 1), and look for the coefficient of the {u — 1)~^ term. After integration, 
this term will lead to the logarithmic divergence. 

With the aid of Mathematica, we find the nested integral can be expanded near the 
horizon as 



u V 



w 



2 + y-4/{l+«,?) (1 



W 



I) 



1 



(2-7;2)2(l- 
2 + m-Hn 



dv 



+ 0{u-l) 



(5.106) 



where Hn{a) is the Harmonic Number defined in (4.78). Using this expansion in Y2, we 
find that near the horizon, 

Y2p{u ^ 1) « 



^ Here the assumption is that C2 > 0, and that gxx{r — > 00) ~ r" with n > These assumptions 
hold for the Schwarzschild AdS metric for any positive p, and for the Dp-Brane metric provided p < 7. 



2/0 I , . / du 



(47rr) 



4(1 - u) 



^ -wt{2 + nT-H„ 

1 — wf 



l + w( 
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(5.107) 



Going back to the general solution for Y2, one thus finds 



log[l-u] <y2b + 



yo 



l-3wl 2 



1-wl 



-wi(2 + i'K-Hn 



(47rr)2 

The requirement of regularity at n = 1 thus gives 



1 + wf 



y2b 



yo 



w{[2 + m - Hn 



1 - Su-fl 
l-wl 



For future convenience, we make use of the identity 



Hn{a) = Hn{a + l)- 



l + a 



to write 



Hr, 



Hr, 



2wl \ 1 - 2wl - ?,w\ 



Substituting this into the equation for ?/26 we find 

2wl 



y2b 



yo 



wi [2 + m - Hn 



+ 0{1). (5.108) 
(5.109) 

(5.110) 
(5.111) 

(5.112) 



Finally, we must apply the Dirichlet boundary condition at u = 0. We proceed as above 
by expanding the integrand of Y2p near u = 0. Mathematica gives 



w 



2 + „-4/(l+«;?) (1 _ ^2 _ ^2) 



dv « {2 + m) + 0{u), 



(2 - i;2)2(l - ^;2) 



so that 



Y2{u ^ 0) «i (4^26 + 0{u)) 



yo 

(47rr)2 



2(2 + m) J 



+ C»(u-2) 



Doing the integral, one finally has 

Y2{u ^ 0) ^ Ay2b - 



yo 



[xQ + Awf (2 + i7r)] +0{i 



(5.113) 

du 
(5.114) 

(5.115) 



and applying Y2{u ^ 0) = gives 



yo 



[xo + 'iwl (2 + 



ITT} 



Using (5.112), 



w\Hn 



2wl 
l + wf 



+ ^(l-3u;?) 



and (5.92) allows us to solve for ^113. The final result is 



W3 



2wi 

1 + 
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(5.116) 



(5.117) 



(5.118) 



This is the main result of this section. To summarize, we have computed the coefficient 
of the term in the sound mode hydrodynamic dispersion relation for a specific class of 



metrics (cf. 5.66). This class of metrics contains two constants which we denote ci and 
C2 . Our expression for W3 should necessarily depend on these constants, but we have 
eliminated ci in favor of T using (5.67) and C2 in favor of wf = Vg due to the relation 

We are now able to plot the sound mode transport coefficients Vg, C; well as this 
new quantity as a function of the parameter C2 to examine how they all behave as 
one deviates from the conformal limit. The conformal limit is the case of C2 = (p + 1)/2, 
giving the Schwarzschild AdS metric. Note that if C2 > (p + l)/2, the bulk viscosity 
becomes negative. Also note that if C2 < p/2, the speed of sound becomes imaginary. 



We only plot the physical region p/2 < C2 < (p + l)/2 in Fig. 5.1 



5.6 Transport coefficients in Israel- Stewart theory 

Now that we have computed the dispersion relation to 0{q'^) in the sound mode, we are 
in a position to examine the implications for both first and second order hydrodynamic 



transport coefficients. As mentioned in Sec. (3.1.3), one of the dominant formulations 
of second order hydrodynamics is the Israel-Stewart formalism |8H [82] . Recently, the 
formulation of second order hydrodynamics presented in [90] has gained popularity, 
though at present it is only applicable to conformal theories. (It should be noted that 
some progress has been made in generalizing the work of [90 to non-conformal theories 
jl28] ). The metrics we consider are not necessarily conformal, and thus we will use the 
Israel-Stewart formulation. 
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The sound mode dispersion relation in Israel-Stewart theory (in the case of met- 



rics with no conserved charge) was given previously in (3.30), and included two new 



transport coefficients: and rn. Comparing (3.30) to our main result (5.118), and 



ehminating t] and ( from the relations (5.85), (5.87) gives the relation 

PVs) 



'1 



1 



(p-1) 



-rji 



p 



{87rT)vUp-l) 



l+vi 



l+2Hn 



l + vi 



0. (5.119) 



As expected, the coefficients t^^ and rn cannot in general be determined separately using 
this method. Still, if one of these coefficients is known, the above relation allows us to 
determine the other. 

Let us now explicitly check that our results agree with other calculations for specific 
backgrounds. The Schwarzschild AdS black hole metric in the near horizon limit was 



given in (4.80). For this metric, the parameter C2 = {p + l)/2. In this case, (5.86), 



( |5.85[ ), ( |5.87| ), ( |5.119D give 

SAdS 



SAdS 
SAdS 



l/4vr, 
0, 

1 



AttT 



P + 1 



P+l 



(5.120) 
(5.121) 
(5.122) 
(5.123) 



Because this metric is dual to a conformal field theory at finite temperature, the results 
for Vs and the vanishing bulk viscosity are exactly in agreement with our expectations. 
The result for r,r is in ag reement with |l08], |IO9]0 This is a non-trivial check on our 
calculation; the cited results were arrived at by completely different methods than those 



we employ here. Furthermore, it should be noted that (5.123) confirms a conjecture 
made by Natsuume |105j ^ 

Finally, we can also consider the case of the Dp-Branes. In the Einstein frame, the 
metric was given in (4.84). For this metric, the parameter C2 = (7p — p^)/(9 — p). Again, 



using (|5.86|), (|5.85|), (|5.87|), and (|5.119), the results for the transport coefficients are: 

(5.124) 

-P 



.DP 



^ In comparing with the results of [109| ■ one needs to employ the identity (5.1101 to see the agreement. 

^ To be precise, the conjecture is confirmed for the case of p > 2; the case of p = 1 should probably 
be checked separately as the derivation of the gauge invariant equations in [2], |101) rely on at least 2 
spatial dimensions. See [129| . where first order hydrodynamics is examined for p = 1. 
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(5.125) 
(5.126) 

(5.127) 



vrT (5 — p) 

The results for v^^ and C^^/v agree with the calculation of |1L)1| . Furthermore, the 
relation between and rn agrees with previous computations for p = 1 and p = 4 in 
[89]. 

This completes our analysis of transport coefficients for this special class of metrics. 



In the process, we have generated formulas for the the speed of sound (5.86), bulk viscos- 



ity (5.87), and a relation between two second order transport coefficients (5.119). These 



formulas are applicable to metrics which obey (4.12), and (5.66), and are generalizations 
of formulas given in jlOlj . 



5.7 Comparison with the shear mode 



Though we have exclusively dealt with the sound mode in this chapter, we now compare 
our results to similar ones from the shear mode. Recall the shear dispersion relation 



(4.1), 



W^(g)shear = -iDr,q^ - iD'^Tshes^rQ^ + 0{q^ 



where 



Ts 



1 



AttT 



(5.128) 



(5.129) 



A formula for Tshear was given in the previous chapter (4.75) and in the paper [T]. This 
formula is applicable to a wide variety of metrics, including the special metrics we have 



considered in this chapter. Using the special metric (5.62 )-(5.66), and the relationship 



(5.67) in this formula yields the result (4.77): 

1 



''"shear 



AttT 



P_ 
C2 



1 



47rT 



Hn 



2wi 



(5.130) 



Before the work of [90 , it was thought that Tghcar = t-k, but as explained in section 
(3.1.3), the authors of the aforementioned paper showed that the coefficient of the q'^ 



term contains not only t^^, but also contributions from (currently unformulated) third 
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order hydrodynamics. We can now determine these unknown contributions for the 
special class of metrics we are considering. 



As in Sec. 3.1.3 let us parametrize 

Tshear = Tyr + A (5.131) 

where A denotes the unknown contributions from third order hydrodynamics. Combin- 



ing (5.130) and (5.119) allows one to solve for A. Evidently, 



The correction A is not universal in the sense that the first term is not present in the 
case of a conformal theory. In the conformal case, A is still not universal, because it 
depends on the number of dimensions of the theory (because Vg depends on p) . In the 
future, when the particular transport coefficients which comprise A are known, it will 
be interesting to see whether there is any universal relationship between these unknown 
coefficients, t-j^, and rn. 



We can easily check that the formula (5.132) reproduces the results in the well 
known SAdS^ metric, which is dual to = 4 supersymmetric Yang-Mills theory at 
finite temperature. In this case, p = 3 and = 1/3. Immediately, we have 



A^^^^.^ = ~. (5.133) 



In [Ml EOj, it was found thalQ 



r^^^^ = ^^M^. (5.135) 
2ttT ^ ' 

It is clear that the relaxation time computed from the shear mode and the sound mode 



differ by the amount predicted by the formula (5.132). 



* Even though the formahsm of Baier et al. [S^ is different than Israel-Stewart, one can check 
that the sound mode dispersion relations coincide in the limit of conformal theories 0). In an 

unfortunate clash of notations, the relaxation time introduced by Baier et al. is denoted by rn- Its 
Israel- Stewart counterpart is r^. 
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5.8 Summary 



In this chapter, we first presented a set of gauge invariant equations for sound mode 



perturbations on a generic black brane type background. The equations (5.42) and 



(5.471 are the main results which pertain to this topic. These equations can be used to 



determine the speed of sound and bulk viscosity for any metric which can be generated 
by a set of minimally coupled scalar fields. In order to determine the dispersion rela- 
tion, one must solve these equations perturbatively in q, applying the incoming wave 
boundary condition at the horizon, and Dirichlet boundary condition at r = oo. 

The gauge invariant equations are quite complicated, and so far a general analytic 
solution has eluded us, though we did present a solution for a particular class of metrics. 
These metrics are not necessarily conformal, and contain an arbitrary number of spatial 



dimensions p with p > 2. Metrics which obey (5.62) and (5.66) have speed of sound 



and bulk viscosity given by (5.86) and (5.87). These results are a generalization of 



the results of |101j . and include both Dp-brane, and Schwarzschild AdS black hole 
metrics. In addition, we have computed the coefficient of the term in the sound 



mode dispersion relation (5.118) for such a class of metrics. 



Information about second order transport coefficients was presented within the con- 



text of the Israel-Stewart theory. In general, a relationship (5.119) between two trans- 



port coefficients r^r and tjj can be determined. In the conformal case of the Schwarzschild 
AdS metric, the relation mentioned above allows the determination of the quantity t^, 
since the coefficient of tjj vanishes in this case. We have verified that our results agree 
with those calculated from different methods. 

Finally, by comparing the sound mode dispersion relation to the shear mode dis- 
cussed in Chapter |4j we were able to determine the contribution of third order hydrody- 
namics to the shear mode (it was pointed out that such contributions would be present 
in [90]). 
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Figure 5.1: Plots of the dimensionless quantities: the speed of sound, ratio of bulk to 
shear viscosities, and u;3(47rr)^, as a function of the free parameter C2 appearing in the 
metric. Each plot corresponds to a different number of spatial dimensions p. Only a 
physical region where both v1 and C, are positive is shown. Note that all of the specified 
quantities vary with C2 and p, and are thus not universal in the same sense as rj/ s. 



Chapter 6 



Conclusion and Discussion 



In this thesis, we started by providing a brief review of the hterature on gauge/gravity 
duahty and its phenomenological apphcations to strongly coupled gauge theories. The 
central focus throughout has been hydrodynamic dispersion relations for perturbations 
of a strongly coupled plasma. In Chapter [3j we gave a pedagogical overview of the theory 
of hydrodynamics and also explained in detail the gauge/gravity duality methods for 
computing such hydrodynamic dispersion relations. In Chapters |4] and [5] we applied 
these methods to a general gravitational dual which is supported by one or more scalar 
fields. 

In Chapter [4] the shear mode dispersion relation was computed up to order (sub- 
leading hydrodynamic order). The equation for Tshear (4.75) is the main result of this 
chapter. This result is in agreement with the original publication [1], though in this 
thesis completely different methods were used to derive it. This is another instance 
where the results from and AdS/CFT calculation agree with a calculation which uses 
the black hole membrane paradigm. 

In Chapter |5] the sound mode perturbations were analyzed in the same type of 
gravitational dual. First, we derived the relevant gauge invariant equations (5.42), (5.47) 
for the perturbations which must be solved in order to compute the dispersion relation. 
These equations are general, and can be applied to any dual metric which is supported 
by (r dependent) scalar fields. Secondly, we solved these equations analytically for 
a special class of metrics, which allowed us to compute the sound mode dispersion 



relation up to order q for such metrics (5.81 ), (5.82 ), (5. 118). Examination of these 
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results within the context of Israel-Stewart theory of hydrodynamics leads to formulas 



for the speed of sound (5.86), bulk viscosity (5.87), and a relationship between two 
relaxation times (5.119). Finally, by comparing the shear mode results of Chapter |4]to 
the sound mode results of Chapter [5] we were able to gain information about the third 
order hydrodynamic contributions to the shear mode dispersion relation (5.132). All of 
these formulas agree with previously known results in the literature, but also offer new 
information about non-conformal theories. 

As mentioned in the introduction, it is desirable to find universal relations among 
transport coefficients, as such relations have phenomenological implications for the 
quark-gluon plasma. Some main results of this thesis such as the shear/sound mode 



dispersion relation results (4.75), (5.86), (5.87), (5.118), (5.119) and the third order 



hydrodynamic contributions to the shear mode (5.132) are applicable to certain non- 
conformal gravity duals in an arbitrary number of dimensions. However, despite the 
fact that these relations are applicable to many theories, they are not universal like r//s 
and the relation presented in |113] . For example, the number of spatial dimensions p 
enters explicitly into our formulas, whereas a universal relation should not depend on 
this quantity. 

Furthermore, we stress that the class of theories examined herein is limited; it seems 
likely that any generalization will explicitly contain the bulk viscosity Q. We were able 
to eliminate this quantity because our analytical results apply only to theories which 



satisfy the algebraic relationship (5.87) between Vg and C,. 

There are many possible extensions of the work presented in this thesis. First, 
the sound mode gauge invariant equations equations can be solved numerically for a 
specified metric and set of scalar field profiles. These equations may be useful for 
phenomenologically based models of the quark-gluon plasma such as those of |43| I116| 
I130j . One chief application of these gauge invariant equations is the determination 
of the bulk viscosity of the dual field theory. Other methods for determining the bulk 
viscosity do exist, in particular the work of |125j is popular in this regard. However, this 
method is currently only applicable to gravity duals supported by a single scalar field, 
while the methods documented in this thesis can be applied to models with multiple 
scalar fields. For example, consider the model of |116] which uses two scalar fields to 
dynamically generate a soft-wall setup. Once a finite temperature generalization of this 
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setup is created, one could use the methods detailed in this thesis to compute the bulk 
viscosity of this phenomenological model. 

In this thesis, we focused exclusively on the hydrodynamic approximation of the 
dispersion relation w^q), which is an expansion in powers of q/T <^ 1. However, it is 
also possible to compute the full dispersion relation w{q) numerically. This was done 
for a conformal theory in |100] . but the analysis could be extended to the special non- 
conformal theories considered in this thesis. Numerical computation of the dispersion 
relation allows one to examine the perturbation's group velocity, (for example, see |13H 
I132j ) . One should check that the theory obeys causality for all values of the deformation 
parameter C2. 

It would also be desirable to increase the generality of the gauge invariant equations 
by including other types of matter such as gauge fields, as this would increase their 
utility. It is a simple matter to add additional matter fields to the linearized Einstein 
equations presented in Chapter |3] since we left the energy-momentum tensor arbitrary 
at this stage and only later specified to the case of scalar fields. 

Another issue which should be investigated is whether there are other special cases 
or regimes in which the gauge invariant equations can be solved analytically. Recently, 
in |133] , the authors find that in the high temperature regime of a single scalar gravity 
dual, their numerical results for the bulk viscosity agree with a simple analytic formula 
to a remarkable degree of precision. This observation may be an indication that a 
general analytical formula for the bulk viscosity exists for single scalar gravity duals. If 
this is the case, one should be able to determine it from the gauge invariant equations 
presented in this thesis. The establishment of such a formula (which was, indeed, 
a primary motivation for the investigation of the work in this thesis) is still highly 
desirable. 

In this thesis we examined the dispersion relations within the context of Israel- 
Stewart hydrodynamics to gain information about the transport coefficients in this 
theory. It would be beneficial to repeat this analysis for the formulation of second order 
hydrodynamics by Baier et al. [90j. We have checked that for conformal theories, the 
sound mode dispersion relations in both formulations of hydrodynamics are identical, 
but it is possible that differences may exist for non-conformal theories. At present, the 
formulation of second order hydrodynamics presented in \90\ have not been generalized 
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to non- confer mal theories, though a first step was taken by Kanitscheider et al. |128j 
(see also |131j ): it would be interesting to combine these results with those presented 
here to see if any new information emerges regarding second order transport coefficients. 
These are all issues on which we hope to report in the future. 
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Appendix A 



General Relativity Reference 



This appendix contains all information necessary to construct the linearized Einstein 
equations (to first order in the metric perturbation) for a general energy momentum 
tensor. The general relativistic notations used herein are those of Weinberg |76j . 

A.l Metric Perturbations 

This appendix is devoted to the expressions for the relevant general relativistic tensors 
when the background metric gf^u is perturbed 

gtj.u = 9^lu + h^y. (A.l) 
Before proceeding one must examine what happens to the inverse metric 

What is Sg^^l By definition, the inverse metric is defined such that 

g^..^" = (A.3) 

{9,. + K,){g^'' + 5g^n = 5^. (A.4) 

Expanding to linear order in h, we have 

g^^ug^' + Kug^p + g^ui^g^ = 5^, (A.5) 
Kug^' + g^,u{5gn = 0. (A.6) 
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Contracting both sides by g'^'^ yields the relation 
Thus 



(A.7) 
(A.8) 



A. 2 ChristofFel Symbols 

The ChristofFel connection is defined as 

1 



r^,. = ij^g^" [df^gau + dyg^^ - d^g^y] . 



(A.9) 



These can be evaluated explicitly for the black brane metric (2.23). The non- vanishing 
results are 

1 



yt 
^ tr 



^ ir 



yr 



-pr 

tt 



^L[gtt{r)], 
T^L[gxx{r)], 

T^dgrrir)], 

-\g''''{r)g[t{r), 
\g"{r)g',Ar)- 



(A.IO) 
(A.ll) 
(A.12) 
(A.13) 
(A.14) 



Here i denotes the spatial coordinates xi...Xp, and we are using the notation Dl to 



denote the logarithmic derivative as defined in the text (3.58). 



Introducing perturbations over the background, one finds the change in the Christof- 
fel symbols: 

<^r^, = Jff^" [d^K, + d,K^ - d^h^,] - h^'' [d^g„, + d,g,^ - d^g^,] . (A.15) 



A. 3 Covariant derivative 



The definition of the covariant derivative changes with the object on which it acts. For 
simplicity, we will act on a vector. Generalization to higher rank tensors is straight- 
forward and is found in any text on general relativity (cf. fTWJ. As in the text, the 



102 

superscript {k) denotes quantities which are kth order in the perturbation. For the 
background metric: 

Vf)?. = 5^^.-r^,6. (A.16) 

Supposing now that we want to examine this quantity under perturbations we straight- 
forwardly have 

'^(V^f^e.) = d^{6^,) - (^r^J^A. - r'^^m). (A.17) 
One application of these formulas in the text is the evaluation of 

6{a('U) = S (5^^V(0)V(0)</.) . (A.18) 
Using the above formula with ^ = du4>, one has 

- g^^-'{6rpdx(P-g^^-'r^^M6^). (A.19) 
A. 4 Ricci Tensor and Ricci Scalar 

By definition, the background Ricci tensor is given by 

4°^ = duTi^ - dxT^^, + Tl^Tt^ - Tl^Tl^. (A.20) 

It's components can be explicitly evaluated for the 'black brane' type metric given in 
(2.23). The results are 



Ft{r) = 5^40) = ^^a, , (A.21) 

F,{r) = /40) = -^a,(^/^Pib.x]), (A.22) 
V 9 

Fr{r) = g^^R(fi} = -^dr {gtt9'"DL[9tt?) + j^dr {g..9''VL[9..f) . (A.23) 

where we have defined the shorthand notation Ft^ Fx and Fr for convenience. 
Combinations which appear frequently in the text are 

Ft-Fx = -^dr {^g9"VL[f]) , (A.24) 
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and 



Ft-Fr 



P 



(A.25) 



where / is defined as in the text ( 3.57| ). 

To first order in the perturbation, the Ricci tensor is expanded as expected 

4i) = d,5T\, - d^5Tl^ + 5Tl^Tl^ + Vl^SVl^ - 5TI,T\^ - TIJT\^. (A.26) 



The Ricci Scalar is defined as 
and in the black brane background takes the form 

Ri.^) = Ft+pF, + Fr 

To first order in the perturbation, 

^(1) = gi^^Ri}) -h^'^R'^!^}. 



"111/ ■ 



(A.27) 



(A.28) 



(A.29) 



One must take great care when raising and lowering indices here. For example, consider 
the quantity R'f^iy This notation means 6 (^R'^q-^^, as follows 



^(1) 



(A.30) 



Notice that R^"^^ / R^^lg'^'^g'"'^ , there are extra terms that come from the perturbation 



metric which raises the indices. 



A. 5 Einstein Equations 

The background Einstein equations are 

4°^ - \R^'^gt^u = -87rGp+2rW. (A.31) 

Here, Gp+2 denotes gravitational constant in p + 2 dimensions, and t'^^J is the back- 
ground energy momentum tensor containing information about what matter supports 
the metric. The perturbed equations are straightforwardly given as 

R^} - \r^^'^9^iu - V = -87rGp+2T«. (A.32) 
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To arrive at the linearized Einstein equations given in the text, examine the quantity 

G'^.^i) = -87rGp+2r(;«, (A.33) 

which means 



-ifj.X 
(0) 



(A.34) 



a.p ap 



-87rGp+2 [Tll^duX + Tl^Q)K\j ■ 

(A.35) 



In deriving this equation, we have used the fact that the Kronecker delta is unchanged 
under any perturbation. Explicit evaluation of the components of this tensor equation 



for the black brane metric (2.23) leads to the equations given in the text (3.50 - 3.52) 



and (3.59 - 3.65) 



We have left the background energy momentum tensor arbitrary at this phase. It 
needs to be computed separately for the particular set of matter fields under consider- 
ation. This is done for scalar field theory in Chapter |4j 



Appendix B 



Techniques used in solving 
differential equations 



This appendix contains some details regarding the techniques used to solve the differ- 
ential equations in the text. 

B.l Reduction of Order 

The technique of reduction of order can be used to help solve second order linear differ- 
ential equations. This technique is useful if one solution of a homogeneous equation is 
known, and one is searching for the second linearly independent solution. 
Consider the homogeneous equation 



and assume that we know one solution yi{x) which satisfies this equation. To find a 
second solution, make the ansatz 



Substituting this into the above equation results in 

v"{x)yi{x) + 2v'{x)y'i{x) + y'({x)v{x) 

+ p{x) [y[{x)v{x) + yi{x)v'{x)] + q{x)yi{x)v{x) = 0. (B.3) 



y"{x) +p{x)y'{x) + q{x)y{x) = 0, 



(B.l) 



y2{x) = yi{x)v{x). 



(B.2) 
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or 

v"{x)yi{x) + v' (x) [2y[{x) + p{x)yi{x)] 

+ v{x) [2/1 (x) +p{x)y[{x) + q{x)yi{x)] = 0. (B.4) 

Using the fact that yi{x) satisfies the difi'erential equation, the last term vanishes, and 
we are left with the first order diflFerential equation for v'{x): 



2/1 W 



= 0. (B.5) 



v"{x) + v'{x) 
This can then be easily integrated to find v'{x) 

v'ix) = i;oexp|-y"(^2^+p(x)j dxj (B.6) 

'"'(x) = — T^expi— p(x)dx> (B.7) 

yi{xr I J J 

where vo is a arbitrary constant. v{x) is then trivially found by performing another 
integration. 

B.2 Variation of parameters 

One fundamental technique used in obtaining general solutions to inhomogeneous dif- 
ferential equations is the method of 'Variation of Parameters'. The method can be used 
to find a general solution to a second order linear differential equation if the solution to 
the associated homogeneous equation is known. 

The theory behind the method can be found in any textbook on differential equa- 
tions. Here, we simply present the essential formulas. Consider a differential equation 

y"{x) + p{x)y'{x) + qix)yix) = g{x), (B.8) 

and assume that the functions yi{x) and y2{x) are linearly independent, and satisfy the 
associated homogeneous equation. That is, 

y'/(x) + pix)y[{x) + qix)y[ix) = 0, (B.9) 

and similarly for 2/2 (x). It can be shown that the function yp{x) is a solution to the 
inhomogeneous equation, where 

/ \ / \ / , / N f y2{x)g{x) 

yp{x) = y2{x) J ^^^^ dx - yiix) J ^^^^ . (B.IO) 



Here, T^(x) is the Wronskian 

W{x) = yi{x)y2{x) - y[{x)y2{x). 



To arrive at the form of the function yp used in the text (5.99), let us define 

2/2(2;) 



Then, 



and 



yp{x) = yi{x) 

One can see that this is equivalent to 

yp{x) = yi{x) J h'{x) j 
by performing an integration by parts. 



h{x) _ 

yi{x) 

W{x) = h'{x)yl{x), 



yi{x)h'{x) J h'{x)yi{x) 
9{z] 



yi{z)h'{z) 



dz dx 



107 
(B.ll) 

(B.12) 

(B.13) 
(B.14) 

(B.15) 



